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The structure of simple decisions is considered in terms of a model 
which composes such decisions from hypothetical elementary decisions. 
It is argued that reaction-time data can be treated by the use of the 
Laplace transform so as to overcome difficulties which negated earlier 
attempts to analyze choice reactions. The general model leads to come 
plex problems which are formulated but not solved. Two special cases 
of the model are worked out, and the statistical problem of evaluating 
the fit of the model is discussed. It is shown that treating decision 
processing as time-discrete leaves the essential features of the analysis 
unchanged. Two experimental proposals, to provide data which should 
be considered in further work on the model, are made. 


I, Introduction. In this paper we propose a model for the way hu- 
man beings organize the decisions required by simple choice situa- 
tions into a collection of component decisions. It is our thesis 
that such an organization of decisions must be reflected in the 
distribution of reaction times and, therefore, that it may be pos- 
sible to infer the organization from the reaction-time distribution. 
Although our thinking derives from empirical studies, we must 
describe this proposal as speculative, for the model is not firmly 
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based on such studies. However, the development of the model 
has led us to suggest two experiments which we believe may help 
to determine what merit it has. These experiments will also help 
to decide whether it is desirable to pursue further work in an at- 
tempt to modify the model to accord better with reality, for we have 
little hope that the particular details of the present model have any 
lasting value. 


Il. Reaction Times. Suppose that a subject receives a stimulus of 
a fixed type at time 0 and responds at time ¢ with a fixed type of 
response. The time interval, ¢, between the stimulus and the re- 
sponse is called the simple reaction time. If the subject is pre- 
sented with one of a set of stimuli and a choice of response con- 
tingent on the stimulus is required the corresponding time interval 
is known as the disjunctive reaction time. In either case, it is 
clear that to obtain stable and readily analyzable time distribu- 
tions it is necessary that the stimulus be simple enough so that 
the mean reaction time is no more than a second or two. Otherwise 
unwanted stimuli may intervene between the test stimulus and the 
response, and the interaction among the stimuli will cause a dis- 
tortion of the time distribution which will be very difficult to 
analyze. 

The study of reaction times, including disjunctive reaction times, 
has a long history in the literature of psychology (cf. Woodworth, 
1938, chap. xiv). In recent years, however, relatively Little in- 
terest has been evident in reaction-time studies. We may attribute 
this loss of interest to two related causes. First, there has been 
a failure to separate the time to make a decision (decision latency”) 
from the other time lags involved in the total process. One attempt 
to make this separation involved measuring the subject’s response 
to a stimulus when no decision was to be made and subtracting 
this time from the time required to respond to the same stimulus 
with the same motor action when a decision was involved. This 
technique has been considered unsatisfactory for the following 
reason. If the subject has no decision to make he is able to bring 
his motor readiness for the specified response to a much higher 


*We use reaction time when referring to the time of a process timed 
from stimulus presentation to motor response; latency when referring to 
times of distinguished parts of such a process. 
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pitch than he can when he is required to make a disjunctive re- 
action; thus, the base time—the time to react in a choice situation 
excluding the time for the decision itself—cannot be equated to 
any simple reaction time. We may conclude that the base time will 
be determined, if at all, only from measurements taken when the 
subject is required to make a decision. 

Second, suppose that in one way or another the pure decision 
latency distribution has been obtained—then what? It is true that 
if these distributions were found to be extremely simple, in that 
they could be well approximated by some class of elementary 
mathematical functions, the separation of non-choice latencies 
(base times) from decision latencies might be an end in itself. If, 
however, the resulting decision latency distribution were of a 
complex character, the challenge to account for it in more primitive 
terms would remain. 

We describe these as related difficulties, for it is not unreason- 
able to suppose that the method used to tease out the non-choice 
latencies (base times) can also be used, or adapted, to decompose 
the decision latencies into more primitive terms. Such a decom- 
position of the observed reaction-time distribution may be an en- 
tirely formal mathematical process with no empirical correlate or it 
may be based on a model which purports to describe the way a 
human being composes the finally observed decision from certain 
more elementary ones. It is with such a model that we are concerned. 

At the heart of our proposal is the idea that the mathematical 
technique of the Laplace transform may be employed usefully in 
the study of reaction times. Since it is unlikely that every one of 
our readers will be familiar with the Laplace transform, we have 
devoted the next section to its definition and to a list of those of 
its elementary properties which we shall need. 


Ill. The Laplace Transform. Let F be a real-valued function of a 
real variable ¢ such that F(¢) = 0 for ¢<0. The real-valued func- 
tion L(F) of the real variable s defined by the equation 


Hea | e~*F(t)dt (1) 
0 


is called the Laplace transform of F. There is essentially no loss 
of information about F in making this transformation [see equation 
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(4)], but because of some of the special properties of the transform 
there is sometimes a distinct advantage to working with trans- 
formed functions. We shall list a few of the elementary properties 
of the transform which we shall need later; no proofs will be given 
for they are well known (cf. Churchill, 1944). 


i. 1} { F (1) Fy (t-1) ar} a L(F LEY: (2) 
0 
it. L (= \- sL(F) + F(0). (3) 
dt 
iii. If L(F) = L(G), then F = G +N, where N is some (4) 


function with the property ie N(é)dt = 0 for all T> 0. If it is known 
that F and G are continous, the N is continuous and so N = 0,i.e., 
Fe=G, 
iv. If a@and 6 are constants, 
L(aF + 6G) = aL(F) + L(G). (5) 


v. If F(t) = Ae~™, where A is a constant, then 


L(F) = ; (6) 
xt! 
IV. The Model. Our proposal is based on assumptions which are 
intuitively acceptable, but which at the moment do not appear to be 
susceptible of direct verification. It is our impression that any 
empirical verification of the model must deal with the full set of 
assumptions rather than with each in isolation. 

Assumption I, It is possible, for a given experimental situation, 
to divide the observed reaction time ¢ into two latency components 
t, and ¢,, called base time and choice time respectively, such that: 

Niet = b+t. 

2. The value of ¢, depends only on the mode of stimulus pres- 
entation and on the motor actions required of the subject. Specif- 
ically, it is not directly dependent on the character of the choice 
demanded. 

3. The value of ¢, depends only on the choice demanded. Spe- 
cifically, it is not directly dependent on the mode of stimulus pres- 
entation or the motor actions required. 

Let the distributions of ¢, t,, and ¢, be denoted by f, fy» and f, 
respectively, Since conditions 2 and 3 imply that the two com- 
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ponent latencies are independent for a fixed experimental situa- 
tion, it follows from condition 1 that 


(=f Klople-nar. (7) 
0 


Our second major assumption concerns only the choice latencies 
and requires the distribution f, to be composed from more elemen- 
tary distributions. The basic idea is that the final decision made 
by a person is organized into a set of simpler decisions which are, 
in some appropriate sense, elementary decisions built into him. If 
such a structure exists in human decision making, it is analogous 
to the structure of a decision process in a computing machine, 
which may be thought of as composed from a set of decisions 
which are elementary relative to that machine, i.e., the elementary 
decision capabilities built into the machine by the engineer. The 
actual organization of these elementary decisions to form a more 
complex one is a function both’ of the individual man or machine 
and of the nature of the decision being made. This is true at least 
of the machine, and we shall suppose it is true of human beings. 
In addition, the breakdown of a complex decision is not, in gen- 
eral, restricted to a serial process where one elementary decision 
is followed by another, for in a machine different portions may be 
simultaneously employed on different parts of the problem. There 
seems every reason to suppose this is also true in a human being. 


We shall describe the organization of decisions by a directed 
graph. (The terms oriented graph and network have also been em- 
ployed in the mathematical literature and the term flow diagram is 
used in connection with computer coding.) A directed graph con- 
sists of a finite set of points which are called nodes, with directed 
lines between some pairs of them. Several examples are shown in 
Figure 1, It is possible, in general, for more than one directed 
line to connect two points, both in the sense that we may have two 
or more in the same direction as in Figure 2a, and in the sense 
that there may be lines with opposite directions as in Figure 2b. 
In this paper, when we use the term directed graph, we shall sup- 
pose that neither of these possibilities is allowed, that is, we 
shall suppose that between any pair of nodes there is at most one 
directed line. 
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FIGURE 1 


We shall employ a directed graph to represent the organization of 
decisions in the following way: At each node we shall assume that 
an ‘‘elementary decision’ will take place, the latency distribution 
governing the decision at node ¢ being denoted by 7, The decision 
process is initiated at node 2 when, and only when, decisions have 
been made at each of those nodes j such that there is a directed 
line from j to 7. We may think of the “‘demon”’ at node 7 waiting to 
begin making his decision until he has received the decisions of 
all the demons who precede him in the directed graph. 

For the directed graphs we shall consider, there will be at least 
one node, possibly more, which is the terminal point of no line; 
these will be the decision points which are activated by the ex- 
perimental stimulus at time 0. There will also be at least one 
node, and again possibly more, which initiates no directed line, 


Sear 
i igen 


FIGURE 2 
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and it is only when the decisions at all these nodes have been 
taken that the motor actions, which signal the subject’s response 
to the experimenter, are begun. It is clear that for any individual 
and for any stimulus situation it is possible to find at least one 
directed graph N and elementary latencies f, which compose as 
described above to give f,. For example, let N have but one node 
and let {,=7~. We shall, however, make stringent assumptions 
about N and f, which, in general, exclude this trivial solution. It 
is some of these assumptions which most likely will be abandoned 
or modified if the present model cannot cope with experimental 
data. 

Assumption II, It is possible to find for each stimulus situation, 
o, a set of stimulus situations, S, which all have the same base- 
time distribution, f,, and an elementary decision latency, f,, such 
that: 

1, o is an element of S. 

2, For each choice situation p in S there exists a directed graph 
N, with the properties, 

a. each of the latency distributions at the nodes is the same, 
namely, fas 

b. the decision time at node z is independent of that at node 
i,j #4, 

c. f, is a composition of N, and 7, (as described above). 

3. Among the stimulus situations in S there is one whose di- 
rected graph satisfying conditions II.2 is a single point. 


In less formal terms, we require that there be groups of stimulus 
situations all of which have the same base-time distribution and 
which can be built up according toa directed graph from elemen- 
tary and independent decisions which all have the same latency 
distribution /,. In addition, among the stimulus situations in this 
class we assume that there is one which employs but a single ele- 
mentary decision. The latter assumption can be weakened, if we 
choose, to the assumption that there is one stimulus situation 
whose directed graph we know a priori, but in what follows we 
shall take the stronger form that the graph is a single point. 


V. Comments. The above assumptions comprise the formal struc- 
ture of our model; there are a series of auxiliary comments which 


are necessary. 
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Even if we were able to show that these assumptions can be met 
for certain wide classes of experimental data, but that in so doing 
we obtain elementary decision distributions 7, which are extremely 
complicated, it is doubtful that we should accept the model as an 
adequate description of the decision process. Equally well, if the 
directed graphs required are excessively complex we should reject 
the model. The hope is that it is possible to subdivide the total 
process into a relatively small set of subprocesses which are 
practically identical. But we do not want to be forced to an analy- 
sis in terms of individual neurone firings. It is probable that As- 
sumption II.3 effectively prevents this extremity by requiring the 
existence of a stimulus situation which demands but one elemen- 
tary decision for its response. 

It is also implicit in our thinking, although not a part of the 
formal model, that the sets S of ‘‘similar’’ stimulus situations will 
include as subsets those experimental situations we naturally 
think of as being similar. For example, suppose the subject is 
presented with n points, one of which is colored differently from 
the others and he is required to signal the location of that one. 
We should want to consider as ‘‘similar’’ the set of these situa- 
tions generated as n ranges over the smaller integers. We should 
probably reject the model if they could not be put in the same set 
S, even if by great ingenuity we were able to find other less in- 
tuitively simple sets of situations for which the model held. 

When the model is applied to experimental data we anticipate 
that the case of the directed graph being a single point will be 
identified with the intuitively ‘‘simplest’’ choice situation within 
the set of ‘‘similar’’ ones. 

In some of the following sections we shall make the following 
explicit assumption as to the form of i. 


Lo afro 20, 

0 uns 
where \ is a positive constant. There are two grounds for sup- 
posing this might be an appropriate assumption. First, let us sup- 
pose that when no decision has been reached by time ¢ following 
stimulation at time 0 then the probability that the decision will be 
reached between ¢ and ¢ + At, where At is small, is approximately 
proportional to At, with a constant of proportionality 4. In this 
case, it is not difficult to show that the distribution of decisions 
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is exponential (Christie, 1952a,b). Whether this assumption is 
correct is an empirical problem, but it must be admitted that it has 
the virtue of simplicity. Second, and probably more relevant, it is 
a relatively common observation that as certain decision situations 
are made more and more simple, the observed latency is better and 
better approximated by an exponential distribution slightly dis- 
placed from the origin (Christie, 1952b; Luce, 1953), The main 
error is generally on the rising limb. If this change toward sim- 
plicity is actually toward a directed graph consisting of one point, 
and if our other assumptions hold, then it seems plausible that the 
elementary decision latency is actually exponential but that the 
observed distribution is smeared by the convolution of the base- 
time distribution and the decision-time distribution. 


VI. The Problem. Let S be a set of choice situations which are 
presumed to satisfy the assumptions of the model, ji.e., S is a set 
of the type described in Assumption II, Let #, denote the reaction- 
time distribution associated with a typical member of S. The 
problem is then to find distributions /, and 7, and a set of directed 
graphs N., where o ranges over S, such that each of the triples 
(7,, £,, N,) when composed according to the assumptions of Section 
IV yields the distribution f,. There may, of course, be no, one, or 
many solutions to the problem, but one hopes that by an appropriate 
choice of S there will be exactly one solution. 

It would appear that if the problem is to be solved in any degree 
of generality, it must be attacked somewhat indirectly. It may 
prove appropriate to solve first the following problem: Given a 
continuous distribution f, find the set of all triples (4, f,, N), 
where /, and 7, are continuous, which satisfy the assumptions and 
which compose to form 7. It seems very plausible to suppose that, 
in general, there are many solutions to this problem. However, if 
f and f’ are two distributions associated with choice situations 
from the same set S, then it will be necessary to accept only those 
triples with the same f, and /, present in both cases. Further 
stimulus situations should serve further to restrict the possibilities. 

These problems will not be attacked, let alone solved, in this 
paper; they appear to be of considerable difficulty. We know of 
only one important lead in this direction, but we have not investi- 
gated it. In recent years, electrical engineers have been concerned 
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with the problem of synthesizing in a systematic manner electrical 
networks to have preassigned transfer functions. If we identify the 
given reaction-time distribution with the transfer functions, the 
graph N with the electrical network, and f, with component char- 
acteristic, there is an analogy between the two problems. This is 
probably worth investigation, but it is almost certain that solving 
our problem will prove to be a major research undertaking. 

To some extent the problem we pose may be simplified by using 
some of our assumptions and the Laplace transform. Let f, be the 
observed distribution of reaction times for a given stimulus situa- 
tion o, then by Assumption II we know there exists a set S which 
includes o and another stimulus situation whose directed graph 
consists of one point. Let f, denote the distribution of reaction 
times in the latter case. From Assumption I we may write 


£(2) = i Roti ddr, 
(8) 


AM=] fOLG-oadr. 
1 i b é 7) QT 
Taking the Laplace transform in each case and applying equation 
(2), 
L(f,) = L(A)LAL) » 
L(f,) ES L(A, L(t) ° 


If we divide the first equation by the second in equation (9), we 
obtain 


(9) 


: (10) 


This is a fairly crucial consequence of our assumptions, for it 
is seen that all mention of the base time has been eliminated. It 
is an equation relating the empirical data to f, and N.. 

At this point we should raise an important practical problem. 
Empirically, one does not obtain estimates of the distribution fs 
but rather approximations to the cumulative distribution 


F(t) -[ f(r) dr. 
0 
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(Throughout we shall use small Latin letters to denote distribu- 
tions and the corresponding capitals to denote their cumulatives.) 
Now, while approximations to F may be reasonably accurate, it is 
well known that numerical differentiation of data tends to magnify 
errors and is, therefore, to be avoided. So the question arises 
whether we can translate our results, in particular equation (10), 
into statements about the cumulative distributions. From equation 
(3) we have 


LG) = sL (FY + FO). 


Since we are speaking of empirical data we may assume F (0) = 0, 
and so equation (10) becomes 


LIF.) Lf) 
os a 11 
L(F,) ~ Lf) oe 


Having eliminated 7, from our discussion, the problem of deter- 
mining it remains. Since our division in equation (11) assumes 
f, is the same in the several cases, it will suffice to determine it 
from any one. The simplest, of course, is the case where the graph 
consists of one point, in which case 


L(f,)  L(Fy) 
jE 4 pe wets Oca be 


L(f,) = (12) 


As an example of how equation (12) may be used, suppose 7, is 
exponential with time constant 4. Then by equation (6), 


Mf 
L(t.) cae ’ 


;t+h 


and so equation (12) becomes 
L(f,) = = L(f,) + L(G) « 


If we make the reasonable assumption that (0) =0, then from 
equations (3) and (5) we find 


ie d 
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Assuming that f, is continuous and that 7, has a continuous de- 
rivative, equation (4) implies 


or integrating from 0 to 7, 
ht bee (18) 


Since f, must be determined from empirical data, it is clear from 
equation (13) that considerable data will be necessary to obtain 
accurate estimates of F,. 


VII. Serial Decision Process. An alternative program to solving 
the general problem discussed in Section VI is to discover the 
consequences of certain explicit assumptions about the directed 
graph N and the elementary latency #,. The results of this alterna- 
tive program will, unfortunately, be much weaker than a solution of 
the general problem, but they may have considerable heuristic 
value. We may choose such extra assumptions on intuitive grounds, 
with the hope that they may be relevant for some experimental 
data. We shall examine two cases which are, in a sense, the two 
most extreme forms of the directed graph N. The first, the topic of 
this section, is the general serial case shown in Figure 3a, and 
the second, which will be discussed in Section VIII, is the parallel 
case shown in Figure 3b, 


a Stimulus @———>-@—-@ -------- -@—_® 


b Stimulus 


FIGURE 3 
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It follows immediately from Assumptions I and II.2.b that the 
observed distribution f, of a serial process having n nodes is given 
by 


t tg to 
d= fff NOG =Gstari vi i 


Applying the Laplace transform to equation (14) and using equation 
(2) we have 


L(f) = L(f,)L(f)” ) (15) 
or dividing by the case n = 1, 
L(f.) mec is) 
ead EG Bn ie a, 16 
Lf) (7) LiF.) (16) 


Equation (16) is the explicit form of equation (11) for the serial 
case. Clearly, if we have given numerical data we may determine 
(possibly numerically) f, for each value of n. 

As an example of how this might be done when we know the 
general form of f7,, suppose f, is exponential with the time constant 
d. In that case, equation (16) becomes 


LE ee cae (17) 
LAO, ace (6 se 
tan 
s 
In Figure 4 we have presented plots of EE ae — for small 
ou d 
values of n. A 
A second equation may be obtained by observing that the mean, 


u,(n), of a serial process with n exponential elementary decisions 
is given by 


py(n) = 1,8) + > (18) 
where ,(5) is the mean base time. Thus, 


y(n) — p(t) = 2, (19) 
We may now use equations (17) and (19) to attempt to decide 
whether a given set of data is adequately fit by the assumptions of 
the model, plus the added assumptions of a serial directed graph 
and exponential elementary latencies. There are serious Statisti- 
cal questions as to how this may best be done, but the following 
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Parallel Case --— 
niT( +2) 
~ (Sens) 


Series Case —— 


S555 


FIGURE 4 
ready method may suffice until the statistical problems are formu- 


lated and solved. From the data we compute as a function of 


es 
LiF, ) 
s; this we may assume is in the form of a plot, which we shall call 


oe A. For each (reasonable) value of n and for some value of = y? 


if 
5 eae eee 
(Fp 
L(F, 
We know from equation (17) that this must be equal to L(F, ) 
assumptions are correct and if the correct value of n has been 
chosen. We thus enter plot A at this point and determine the value 
of s. Since we selected \ = 2s, this determines A. But equation 
(19) presents a relation between the observed means, A, and n 
which will be satisfied if our assumptions are valid. We choose 
the value of n such that the error between the observed means [the 
left side of equation (19)] and (n —1)/A is a minimum; this yields 


say — =e find in Figure 4 the corresponding value of 


if our 
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the best possible fit at the point s/\ = 1/2 for the model with the 
added assumptions of a serial graph and exponential f. Using 
these values of A and n, one may add the theoretical curve 


1 <<; VS. 8 to plot A, and a comparison between the two 

S} ~ 
(x + 1) 

curves will give some indication of the adequacy of the assump- 


tions. Clearly, a less subjective criterion of the quality of this 
fit is needed, 


VII. Parallel Decision Process. If we suppose that the n elemen- 
tary decision processes are carried out in parallel (see Figure 3b), 
the choice latency distribution is the distribution of the largest of 
nm selections, one from each of the elementary distributions. This 
is known to be given by 


phage PF 
= At 
att hi AD ; 


which in the case when all the elementary distributions are the 
same, namely F’, , reduces to 


nf, (t)F(2)1"~* . 


If we denote the observed reaction-time distribution for the parallel 
case by g_, then it follows from equation (7) that 


g,(t) = f (rnp (t - 1) F(t - 1) \-*dr . (20) 
0 


Applying the Laplace transform and equation (2), 
L(g.) = Lf, LnfF,*?) « (21) 


As before, we may divide by L(g,) to eliminate L (f,) - 
To proceed further, we assume f, is exponential, then 


Lp ea ey t emte“M[L = eA! ae , 
0 


co nal 
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0 k=0 
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To evaluate the above sum, consider the function 


n—1 


D(z) = (27) (~1)Fa* wt 3 oe (1-2)?! 
x & ; 


Observe that 


1 Fe fa 
af (adde= nD CEN of aaeeee 


n—1l 


= ("E')(- 1)* —_=__ ’ 
5 Lear 


=Li(nfcRike =, 
and that 
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0 0 
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where B(m, n) is the Beta function and I'(n) is the Gamma function. 
From these results we easily obtain 


> 


L(g,) nB(> +1, 
L(9,) " B(} + lf 1) 


nir(E + 2) 
: ee (22) 
ry +n + 1) 
nit (> + 2) 
In Figure 4 we have also presented plotssof — aya 


for small values of n. as 1) 
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The mean of the parallel process can be shown to be given by 


2 Awana 
4 (n) = pw, (0) + = Me - and thus we have, as in the serial case, a 


t=1 
second relation which must be met 
Liceah ih 
n(n) ~n,(l)=5 "5. (23) 
The procedure for curve fitting is the same as described for the 
8 


serial case except that a 1 seems to be a more favorable place 


to enter the graph than is : Soe 


IX. Model Selection. Without a solution to the general problem 
described in Section VI, there arise statistical problems as to how 
well a particular set of assumptions, such as serial directed graph 
and exponential f,, fit the data and whether another set of similar 
assumptions is better or not. In addition, within any one set of 
assumptions there are undetermined constants, such as A and n, 
and there is a question as how best to choose them. We have indi- 
cated one procedure (end of Section VII) to determine the con- 
stants, but it is almost certain that such an ad hoc procedure is 
not optimal. 

The difficulty of making a selection among different sets of as- 
sumptions is evidently quite serious for it can be seen from Figure 
4 that for almost any small value of 7m in one there is an n’ in the 
other such that the two curves are fairly similar. Presumably, any 
other directed graph will produce curves which, in some sense, Lie 
between these two extreme cases. Thus, the shape of the empirical 
data curves will not be extremely revealing of the proper directed 
graph to use—an unfortunate situation. 

It is clear that there are a number of difficult statistical prob- 
lems here, but in all likelihood it will prove to be more efficient 
first to do some experimental exploring using subjective judgments 
as to goodness-of-fit before trying to formulate and to solve the 
statistical problems. 


X. The Perceptual Moment. In Section Il we remarked that in 
reaction-time studies the mean reaction time should be of the order 
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of one second if unwanted interactions with other stimuli are to be 
avoided. This means that the data will be in a range where certain 
peculiar phenomena have been observed. To explain these ob- 
servations, it has been proposed that a subject processes informa- 
tion very rapidly at certain discrete times and that he is in a re- 
fractory period between them. The period from the beginning of 
one such hypothetical event to the beginning of the next has been 
termed the perceptual moment (Stroud, 1949a,b). Unfortunately, 
relatively little direct experimentation has been conducted on this 
problem, and so it is not possible at this time to give a formal 
statement of the properties of the moment. Indeed, there are in- 
vestigators who doubt its existence. In the case that it does exist, 
our analysis will be applied to situations where it most probably 
will have an effect. It is, therefore, of interest whether the analy- 
sis can be adapted to cope with it. In this section we shall make 
a simple hypothesis as to the nature of the moment, not with any 
belief that it is correct, but only to indicate that the general fea- 
tures of the analysis remain unchanged. 

Let us assume the moment is of fixed duration, say 5 seconds, 
and that while a person may receive information at any time during 
that period it will only serve as a stimulus at the end of the pe- 
riod, Furthermore, we will assume that all intermediate (elemen- 
tary) decisions occur at multiples of 5. Since we may assume that 
there is no correlation between the stimulus presentation and the 
timing of the moment, we may assume the stimulus is presented ac- 
cording to a uniform distribution fA in the interval 0 to 6. This 
assumption may be inappropriate, for it may happen that a person 
is only able to assimilate information during part of the moment; 
we shall return to this point later. 

The question now arises as to the discrete form we should as- 
sume for the elementary decision process. In the continuous case 
we took it to be exponential, and so we shall use the discrete 
analogue. We assume that if no decision has been reached by the 
ath moment following the presentation, i.e., at time 76, then the 
probability of a decision in the zth moment is AS. If we call the 
probability of a response by the zth moment P,, then 


P, = P,_,+ U-P, 308 ,; (ou 


= (1-d5)P, , +28. 
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With the initial condition P, =0, the difference equation (24) is 
solved by 
P,=4—(1~)8)'. 


The probability of a decision in the zth moment is obviously 


[bP oe Ade; 
hence, we have 
AO( Leno ee (25) 
as our distribution Bx 
If we replace this discrete distribution, equation (25), by a 
continuous one ®, which has rectangles of width ¢ and height 
rb (1 ~— dS)?! 
€ 
the limit as « —» o this becomes the discrete distribution. 
Let the base-time distribution be denoted by 7, as before, then 
the observed data in the discrete serial case is given by 


f(t) = Jim [ alee fy (t,) Mig - t1) Ot, — ty)... 
0 0-0 
O(t-¢ 


centered about the point 26, then it is clear that in 


)dt,...dt 


n+1 n+1 
Applying the Laplace transform and using equation (2), 


L(f,) = lim L(f,) L(A) L (®,)" = L (f,) L(0)| Jim L @ap Sa!) 


Observe, 
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But, 


so, 
lim L(®,) co 
; Se Sv Fas Gia Z 
an a ere 


Substituting in equation (27) and dividing by the case n =1, we 


have 
n—1 


ify) Speen (28) 
L(f,) -| ee ee etd ee 


which is the crucial equation for the discrete serial case. The 
mean of the discrete distribution /, is given by 


co 3 1 

Y. as(1-r8)Ph =>. (29) 
i=1 

Thus, the relation between observed means is 


n-1 


py (n) ~ 1,11) = =—. (30) 


Now, if we know the value of 6, i.e., the length of the moment, 
then these two sets of equations may be used in exactly the same 
fashion as were equations (17) and (19) of Section VII. We have 
no theoretical value of 5, so it will be necessary to perform in- 
dependent measurements of it. It is clear that if the perceptual 
moment is a real phenomenon it will be important to ascertain its 
properties prior to analyzing experiments on reaction time. 

One further comment of some interest: If we ignore f, and let 
n =1, the convolution of A and ®., when « —> 0, is a step function 
such as that shown in Figure 5, The convolution of this function 
with 7, , for reasonable f, , will serve to smear the steps but it will 
not utterly destroy them. Smearing will also result if n is larger 
than 1, the amount depending on the value of n. Thus, if our as- 
sumption as to the moment is roughly correct, we should expect, 
at least for comparatively simple situations, to find the observed 
latency distribution somewhat lumpy. Indeed, in the literature (cf. 
Woodworth, 1938) it has been remarked not only that the data are 
lumpy but that there is an oscillation superimposed on the distribu- 
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Convolution of h and ¢, 


§ 26 36 46 55 
FIGURE 5 


tion curve. This effect could easily be obtained analytically if we 
were to assume A uniform over only a small portion of the interval 
0 to 6, in other words, if we assume the vast majority of the mo- 
ment is truly a refractory period during which there is no intake of 
information. These considerations bring out even more strongly the 
need for comprehensive experiments to determine the properties of 
the moment. 

We shall not attempt, as before, to study the parallel case. The 
reasons are that the mathematical problem is rather complex and 
with so little information on the nature of the moment it hardly 
seems worthwhile to carry out the analysis. Furthermore, we are 
of the opinion that it is unlikely that information accepted in dif- 
ferent moments is dealt with other than serially. It may happen, 
however, that the information accepted in one moment is processed 
in parallel. The latter remark is a possible hint for developing an 
explanation of the effect of changing the number of ‘‘psychological 
dimensions”’ in an information display. 


XI. Experimental Proposals. The key assumption in our analysis 
is that elementary decision processes can be found of such a sort 
that complex decisions can be built up from them in a way which 
leaves their characteristic 4 value invariant. One should Like to 
present experimental subjects with stimuli which vary in several 
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dimensions but for which decisions on each of the dimensions have 
identical time characteristics. If one uses conceptually different 
dimensions, we run into the difficulty of possibly introducing 
several different \ values. If we use several objects with the 
same dimension relevant for each and with identical characteristics 
in every other respect, we have the difficulty that the reception of 
the stimulus may not be unitary, but broken down into several 
parts separated by receptor orienting acts such as eye movements. 
The first of the two following proposals suffers from the latter 
difficulty; the second from the former. 


1st Experiment: Digit Difference Perception 


Stimuli: White 3" x 5" cards with a triple-spaced typed, horizontal 
row of vertically aligned pairs of digits, 0 and 1, on each. The 
number of pairs per card to vary from one to sixteen. On each card 
either one pair or no pairs will be unlike digits, i.e., (0,1) or (1,0); 
the remainder like pairs, i.e., (1,1) or (0,0). The place of the un- 
like pair in the series of pairs to vary from the initial to the final 
position. Cards with the unlike pair in each of the positions from 
one to n will be included in the set with equal frequency, and 
cards with no unlike pair will be included with the same frequency. 
The assignment of (1,1) or (0,0) to the remaining places will be 
made on an equiprobable random basis, and the choice of (0,1) or 
(1,0) for the unlike pair will be made on the same basis. 


Responses: Experimenter will announce prior to each stimulus 
presentation how many pairs the card to be shown bears. Subject 
will respond yes or no, depending on whether the card does or does 
not bear an unlike pair, by pressing the appropriate one of two 
keys. The subject will be told that an unlike pair in each of the 
possible positions, including in no position, are equally likely 
events, and will be instructed to read the lines of pairs from left 
to right. The data of primary interest will be the latencies of the 
mo response to the cards which bear no unlike pair and the la- 
tencies of the yes response to the cards which bear an unlike pair 
in the nth position. 


Apparatus: 1, Stimulus cards as described above, 
2, Light projector with fast shutter, 
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3. Three telegraph keys: (a) for the subjects to rest 
their fingers on prior to response so that the re- 
sponse will always start from the same situation. 
(b) for yes responses (c) for no responses. 

4, A buzzer of 4% sec duration as a warning signal to 
be sounded ending 1 sec before shutter opens to 
illuminate stimulus. 

5. Recording chronoscope accurate to at least + 10 
millisec. 

6. Timer for ready signal and shutter operation with 
silent starting key for the experimenter. 


2nd Experiment: Multi-attribute Perception 


Stimuli: Ten decks of 32 cards each to be prepared using two 
values on each of five attributes according to the following scheme: 


Attribute Values 
1. Number of spots 2; 3 
2. Color of spots Red; black 
3. Shape of spots Round; square 
4, Arrangement of spots Horizontal line; vertical line 
5. Background color White; green 


Responses: Experimenter will announce what pattern of attributes 
is to be responded to positively prior to each stimulus presenta- 
tion. Subject to make a yes or no response by pressing the appro- 
priate one of two keys as exemplified below: 


Experimenter Says Stimulus Presented Sto Respond 


1, Round red Two black squares in No 
horizontal line on white 
card 
2. Vertical line of squares Three red squares in Yes 
on green card vertical line on green 
card 


The instruction-stimulus pairs which call for a negative response 
should be half of the total number of stimuli presented in each 
attribute-pattern category so that the uncertainty of response prior 
to stimulus presentation will be equalized at the maximum. The 
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data of primary interest will be the latencies of response to the 
set-stimulus pairs calling for a yes response. 


Apparatus: Same as for the first experiment except for the stimulus 
cards. 
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This paper is a continuation of some discussions published previously 
(Bull. Math. Biophysics, 16, 317-347, 1954). A different example of a 
primordial graph is Ae and discussed. Then six different types of 
geometric transformations, similar in some respects to the one discussed 
in the previous paper, are proposed. Another set of seven different 
transformations which emphasize the central role of the circulatory sys- 
tem is described. A theorem about the ‘‘lumping’’ of several points of a 
graph into one point is proved and possible biological aspects of it are 
discussed. 


This paper is a continuation of a series of preceding ones 
(Rashevsky 1954; 1955a,b,c; referred to correspondingly as I, II, 
III, IV) and presupposes familiarity with those. 

In I we have discussed a primordial graph which was used merely 
for purposes of illustration without any claims as to its actually 
describing a primordial organism. However, as we pointed out (I, 
p. 328), some features of the graph are, if not dictated, at least 
suggested by actual biological considerations. As has been re- 
marked in I (p. 328), the determination of the graph of an organism 
is ultimately a matter of experimental study. Partial graphs of 
organisms have been determined experimentally and are well known 
in physiology. We may cite as an example the graph which repre- 
sents the interactions of different endocrine glands or the graph of 
some biochemical reaction as for example the graphs of the glyco- 
lytic cycle and of the ‘‘citric acid cycle’’ shown on pages 25 and 
26 of B. Scheer’s (1948) ‘‘Comparative Physiology.”’ 

If we could determine experimentally the complete graph of any 
organism, we could then use it as a basis for a bio-topological 
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theory, deriving various bio-topological properties of the organism 
from the topology of the graph. Because of the tremendous com- 
plexity of the graphs of even the simplest organisms such a pro- 
cedure is not feasible. Therefore, there is no choice left but to 
postulate a primordial graph and a transformation and see whether 
the conclusions derived from the postulates are in agreement with 
observations. Such a postulational program leads us, as we have 
seen in IV, to different ‘‘systems of abstract biologies.”’ 

Leaving the discussion of the primordial graph still on a purely 
abstract theoretical level, let us‘ nevertheless discuss some possi- 
ble modifications and improvements over the graph proposed in I. 


FIGURE 1 


There is one process in even the simplest organism which oc- 
cupies a somewhat special position both biologically and topologi- 
cally. Even a simple cell is not homogeneous, and the different 
biological and biochemical reactions take place in different parts 
of the cell. A connection between such processes is made possi- 
ble only by transport phenomena. The principal mediator of trans- 
port is diffusion, though sometimes convection may play a part. 

To indicate that a reaction f, cannot take place without the re- 
action 7; previously taking place, but that /, also needs a transport 
f,, of reaction products from the site of f; to that of f,, we symboli- 
cally made converge on the representative point of /, two arrows: 
one connecting /, and /,, the other connecting the representative 
point /,, of the transport to f,, (I, p. 327), as illustrated on Figure 
1. With this method of representation the assumption is made that 
in the primordial graph a biological function takes place only if 
all the immediately preceding biological functions take place. 
This does not need to be the case in the transformed graph where, 
according to theorems 1 and 2 of IV, arrows from points which all 
represent the same biological function in different component pri- 
mordials may converge on the same point. In fact, the biological 
interpretation of the above theorems as given in IV implies that in 
the transformed graph a biological function takes place without 
necessarily al] immediately preceding ones taking place. What is 
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needed is that at least one of each biological functions which pre- 
cedes the given biological function in the primordial should be 
taking place. 

The method of representation illustrated in Figure 1 has obvious 
disadvantages. The direct connection f; — f, does not mean the 
same thing as the direct connection f,— f,, between representa- 
tive points of two reactions of which one follows the other without 
transport. Yet topologically the two connections are of the same 
kind. If transport takes place by convection, then the same con- 
vection current designated by a representative point f, (Figure 2) 


FIGURE 3 


may carry the end products of reaction /, to the site of f,, as well 
as the end product of reaction f, to the site of f,, resulting in a 
diagram shown in Figure 2. That graph does not indicate clearly 
that f,, stands for the transport of two different substances by the 
Same means. 

We shall, therefore, use a different method of representation 
which obviates those shortcomings. The proposed new representa- 
tion will have particular advantages when we discuss the trans- 
port in a differentiated organism represented by a transformed 
graph. 

We represent the general notion of transport by a representative 
point f,. (Figure 3). This transport may carry different substances 
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from different sites and involving different reactions f,. If the end 
products of the reaction /, are transported to the site of f,, we rep- 
resent the transport of that particular product by a point /;,,. 
Since that particular phase of the transport cannot take place with- 
out the existence of the general transport mechanism f,, we there- 
fore make the connection f/f, —> f,,, 2S in Figure 3. Then we make 
the connection /;,, —f,- If there is a transport between two 
other reactions, represented by the points f, and /,, we represent 
it by the point 7, and again make the connections f,, — f,,, 
and is Pees irs 

A situation may arise where a reaction represented by 7, will 
take place only if end products of two other reactions 7, and 7, are 
transported to the site of f,. The transport of the two kinds of end 
products is represented by the point f,,,, (Figure 3), with the 
connections: f, — f,,rwi Ty — he pee 
fav Tw ae hy 

In line with this method of representing transport phenomena, 
let us consider the following possible primordial graph or rather 
the following partial graph of the primordial. For convenience we 
use numbers for the subscripts of the points f; The biological 
function represented by each point is explained in the text. 


— f,. 


Ingestion f, (Figure 4) is followed by digestion f, which leads 
on one hand to defecation f, and on the other to absorption f,. The 
absorbed materials are transported via f,,, to the sites of syn- 
thetic processes f,. Then the synthesis of digestive enzymes, 
represented by /,, follows via transport 7/,.,. Those enzymes are 
transported via f,,., to the site of secretion, represented by f,, and 
digestion f, again follows. 


On the other hand, some of the synthetized products are trans- 
ported via f,,., to the site of the catabolic processes, which are 
represented by 7/,. Products of catabolism are transported via lave 
to the site of elimination of waste products, and that elimination, 
represented by f,, takes place. Catabolic processes result in the 
liberation of energy, represented by fio, Which in turn provides the 
possibility of transport f+ On the other hand, after a transport 
ler11 the catabolic reactions give rise to the production /,, of 
CO,. The latter is transported within the cell via firrigs and this 
eventually results in the elimination of CO,, this process being 
represented by /,,. 
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The intake of O, from the outside, represented by f,3» results in 
a transport of O, to the sites of different reactions involved in 
catabolic processes (Scheer, 1948). Liberation of energy combined 
with anabolic processes as well as with other biological functions 
results in the process of multiplication, which is not indicated in 
Figure 4 to simplify the latter. 


FIGURE 4 


All the above discussed biological functions belong to the class 
of vegetative functions. A partial graph of the animal functions 
may be represented as follows: : 

A stimulus /,, combined with energy release fio Tesults in some 
kind of protoplasmic motility f,, which leads to the ingestion hye 
A different kind of stimulus f,,, combined with energy release f,,, 
results in a different set of motility processes /,, which produce 
an avoidance reaction /,,. 

The primordial graph suggested here is different from the one 
used in I and possibly somewhat more realistic. It is, however, 
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still as schematized as the other. We still omit a tremendous num- 
ber of biochemical reactions of an enzyme nature which are con- 
nected to the different biological functions shown in Figure 4. 
The omission of a number of points, while making the graph much 
simpler, is not necessarily objectionable in principle. What is left 
after such an omission is a partial graph, and if we use the trans- 
formation 7 discussed in I, we can apply theorem 2 of I (p. 338) 
and thus derive some properties of the corresponding partial graph 
of higher organisms. 

A more serious objection is the circumstance, mentioned already 
in I (p. 328), that some points of the suggested primordials stand 
actually for possibly a very large number of points that are inter- 
connected in an intricate manner. Thus point 7, represents actu- 
ally a very complex set of reactions in which very large numbers of 
different organic substances are synthetized. If those different 
points are specialized in different tissues, then the permissibility 
of such a ‘‘lumping’’ of whole partial graphs into one point may be 
seriously questioned. This problem will be discussed below after 
a study of different transformations. 

In discussing the transformations of the primordial graph, .we 
considered in I that all points of the primordial are divided in two 
classes: specializable and residual points. We considered the 
situation in which a biological function when specialized by one 
tissue is lost in all other tissues. While this is true with respect 
to some biological functions, it is definitely not true for all of 
them. In IV we pointed out that such biological functions as trans- 
port, intake of oxygen, elimination of CO,, and elimination of other 
waste products are retained by all cells, although in higher or- 
ganisms specialized organs carry out those biological functions for 
the organism as a whole. Such biological functions we shall call 
specfalized residual. A closer examination reveals that they are 
not aS numerous as the above list suggests. 

While it is true that every cell of a higher organism takes in oxy- 
gen, this intake is different from that of a protozoan because the 
latter receives the oxygen from the inorganic environment directly, 
whereas each cell of a higher organism takes in oxygen from blood 
which, while being exterior to the cell, is nevertheless an integral 
part of the organism. The specialized organs, such as gills or 
lungs, retain the property of taking up oxygen directly from the ex- 
traorganismic inorganic environment. The same holds about elimi- 
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nation of carbon dioxide and of other waste products. A somewhat 
different position is occupied by digestive processes. In the high- 
est animals all digestion occurs extracellularly in a specialized 
organ (the intestinal tract) so that practically all other cells have 
lost the ability of intracellular digestion. In lower forms, such as 
Platyhelminthes, Nemertea, and others, which do possess a more 
or less differentiated and specialized gastrointestinal tract, di- 
gestion takes place partly intracellularly, partly extracellularly. 
Thus in spite of specialization not all cells lose the intracellular 
digestive functions. 

We must, however, keep in mind that even intracellular digestion 
is extra-protoplasmic, occurring in the digestive vacuole. Techni- 
cally the digestive vacuole is inside the cell. Yet it is not part of 
the living protoplasm. In bacteria digestion takes place partly 
extracellularly. Thus we may say that what is actually special- 
ized in a higher organism is the secretion of digestive enzymes, 
represented by /, in Figure 4. The process of digestion is basi- 
cally an external process. All other cells of the organism lose the 
ability of secreting digestive enzymes though some specialize in 
their production. Thus f, is a specializable point, and it may be 
said that f, and f, always specialize in the same organ. 

What actually happens in the case of oxygen intake or carbon 
dioxide elimination is that those two biological functions are spe- 
cialized by an organ consisting of several tissues (IV, p. 213) and 
that first oxygen is carried by the intracellular transport of those 
tissues then taken over by the specialized general transport, if it 
exists, then taken over again by the intracellular transport of each 
cell and carried there to the site of f, (Figure 4) or any other ap- 
propriate site. It is, therefore, more appropriate and rigorous to 
designate f,, as the point representing intake of oxygen from in- 
organic environment, f,, aS representing elimination of carbon 
dioxide into inorganic extraorganismic environment, and so forth. 

The organ par excellence of the specialized transport is the cir- 
culatory system proper. However, the lymphatic system as well as 
the interstitial fluids may, in general, be considered as part of the 
specialized transport. There remains the important question as to 
whether the functions of the vascular system, the lymphatic sys- 
tem, and the interstitial fluids are to be considered as specializa- 
tions of three distinct transport functions present in the primordial 
graph, or whether they are biological functions which correspond to 


120 N. RASHEVSKY 


different subsidiary points in the specialization of a single biologi- 
cal function, that of general transport. The latter view appears to 
be more attractive. However, only further detailed studies of dif- 
ferent possible transformations, some of which will be presently 
discussed, will enable us to decide the above question. In the 
present paper this question is left open. Whenever we shall speak 
in the present paper of specialized transport, we shall explicitly 
understand the cardio-vascular transport. But we do not commit 
ourselves here as to whether to consider the lymphatic system as 
being represented by subsidiary points or not. 

Thus, the only specializable residual points of the graph of 
Figure 4 are the transport points f,, and the /,,.,’s. Their speciali- 
zation may be described in the same manner as that of other points. 
In other words, if for example the sth component primordial graph 
(IV) specializes in all or some of the transport points, it loses 
some of the other specializable points, and a corresponding number 
of subsidiary points is added to the specialized residual points 
which represent the transport. Although biologically apparently 
either all or no transport functions specialize, yet mathematically 
we may consider the case in which only some of the points f,,, 
fz7% Specialize. 

In the zth component primordial graph which specializes in some 
non-transport points only, we now connect each specialized point 
f{? not necessarily directly with those remaining points Se of 
each graph & (k = 1, 2,..., 2-1, ¢+1,..., m) which correspond 
to those points fh? of graph z to which ae is connected. We do 
this only if f“* is not a transport point. If it is, in other words, if 
[ie fai ae or ee feos, where the dots stand for some aie 
letters, and if no transport points are specialized, then, in line 
with what was said above, we connect /(? _ to all fs in such 
a manner that 


1 > 12>. (1 


would be a directed way from fo to foe . Or if See iil then 


we connect /,,, to all ey ry in Rise a manner as to hae. a di- 
rected way 


ies 4 alee ek ee (2) 


If, however, the transport point fir... © f...7, 18 Specialized, say 
in the sth component primordial, then we connect te (resp. 
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fetd) first. to feom (resp. f‘*,,,) and to all its subsidiary points 
s 

fir ...,a (PSP. A p7,0)3 [A = 1,2,... EP 5 or & = 1,2,..., 7677] 
and then connect f{3) (resp. f‘*?,,,) and all the subsidiary points 


to each ees (resp. oe) in such a manner that a set of directed 
ways 


(i) (s) (k) k 
fir... fir... > fir... and ess: OOS eon Ieee (3) 
respectively 

(i) ( (k) ec k 

Cte peel and Pee Patel re (4) 


would be formed. Thus each specialized biological function which 
is preceded or followed by transport is connected to others via 
transport in its own tissue, systemic transport, and transport in 
each of the kth tissues. 

We now are ready to formulate a transformation different from 
transformation T discussed in I and III. Before we do that, how- 
ever, we Shall make a digression and mention six variations of the 
transformation 7, the possibility of which was already mentioned 
in I (p. 332) and which we shall denote by (EES EY AL) menor SEU 

Instead of connecting the subsidiary points in this manner (I, p. 
334; III, p. 112): 


) _, 4@ ; ©. 
fs fh ecty ile (5) 
we may connect them thus: 
(i) (i) Cis 6 
Nien cor ace (6) 


Or we may connect them thus: 


(i) (i) (Whe ” 
if DOO OS ON (7) 
Finally we may not connect the subsidiary points to each other at 


all, connecting them only to those points to which fn is con- 
nected. We shall symbolically denote this by 


Hered. Ae , (8) 


Instead of connecting a specialized point with its subsidiary points 
“in series’’ as indicated by (5), (6), or (7), we may connect them 
‘‘in parallel,’’ and this in three different ways. 
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Either we make the connection 


(i) ed On (i) (2) oer 
fi fie ele ious SNES. : - 


or the connections 


; Dire 
fom Te) lee eye oe 
or, finally, 
Ct) es @MwfO,. . 11 
ie i; pooh ine ( ) 


We now define the following transformations: 


T°?) is identical with T, except that (6) is substituted for (5). 
T“ is identical with T, except that (7) is substituted for (5) 
T‘*) is identical with 7, except that (8) is substituted for (5). 
T° is identical with 7,-except that (9) is substituted for (5). 
T‘®) is identical with 7, except that (10) is substituted for (5). 
T*™ is identical with 7, except that (11) is substituted for (5). 


Now we define a transformation 7 in line with the discussion 
preceding the definitions of T™...7™. 

Let n biological functions /, all arbitrarily numbered from 1 to n, 
become specialized in m groups of n,, n,,...,n,, points each. We 
have as before 


>on, = 03 n>. (12) 

Let group ¢ specialize in n, biological functions f,, where / desig- 

nates one of the n, points. Since, in general, n may include some 

of the transport functions which are both specializable and resid- 

ual, therefore, now the number of points lost in group 7 is not nec- 
essarily n~1n,. We shall denote it by 7,. We have 

n.<n-n,, (13) 


t u 


because transport points are included in the n points, but they are 
never lost. 


Let f, designate one of the 7, lost biological functions.* 


*In I, p. 111, it was stated that the n— n, points f, are divided into 
m- 1 groups, each consisting of Try Moyeee ss 4) Mey yyeeeyM, points. 
This is an error. For the case, discussed in I, the n—n. lost points fie 
are divided into m— 1 groups, each consisting of n— nine Noy eee 
M43 M—N:445 ++-N, points. This relation also holds now, but only if 
none of the transport points are specialized. 
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The transformation rules for 7 are: 


(8) 
rT; 

(8 
rs 


(8) 
rs 


(8) 
ie 


(8) 
Ts 


Draw m distinct primordial graphs. The number m is a param- 
eter of the transformation. 

Remove from graph z (¢=1,2,...,m) those m, points which 
correspond to its lost biological functions. 

Cae if {e is a remaining point of one of the m—-1 graphs 
k (k = 1,2)... 59,4) 449+++ 5m), Which corresponds to a point 
i to which [? is connected, and if fe is not a transport 
point, then connect /{? to each such f‘*), preserving the di- 
rections of connections between ae and f®, 

(2). If an f< is a transport point fie. & f 2. -and if=tio 
transport points are specialized, then connect ews (resp. 
f,,;) to all points /{P (resp. f,,,), preserving the direc- 
tions of connections, in other words so that expression (1) 
respectively (2) would hold. 

(3). If an i is a transport point f{? (resp. f .,) and if 
that point is specialized in the sth graph, then connect Lee 
(resp. f,.,) to f{ (resp. f%,,) and the latter to each 
f® - (resp. boop (k = Sy ne 7 — 4» POP) AEE to 

. t 
have a directed way from f{p___ to f{)__ (resp. from (0 to 
fT) 
Choose a distribution of 7, elements into n, sets, so that the 
lth set has a number wi?) of elements with 


nN. 
t 


>, vO =%, . (14) 
=1 
Some of the yi? may be zero. This will always happen when 
n,>n;, This arbitrarily chosen distribution plays the role of 
another parameter of the transformation. A third parameter is 
the choice of the n,’s. 
In graph < add 7, subsidiary points and divide them into n, 
sets containing each pe points. Connect each of the oss 
points of a set to all other points to which the specialized 
point f{? is connected, in the same way as f<? is connected 
to them, preserving the directions of connections, and also 
make the connections 


1 12, > +120 (15) 
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7) Connect each point ff? (r= 1,2,... ,v$) to a residual graph 
in the same manner as tees is eee tied to its own residual 
graph. 

Application of OSL to the component primordial graph in which 

transport points are specialized insures connections (3) and (4) 


automatically. 
In a similar manner as we obtained transformation 7 ... 7 
from 7, we can obtain six other transformations T® ...7° from 


T°), We list them here: 


T°) is obtained from 7“ by substituting (6) for (15) in TC. 
T“) is obtained from 7“® by substituting (7) for (15) in T°), 
T“ is obtained from 7“ by substituting (8) for (15) in T®), 
T“) is obtained from 7“? by substituting (9) for (15) in Te), 
T“®) is obtained from 7‘® by substituting (10) for (15) in Ts, 
T“*® is obtained from 7‘® by substituting (11) for (15) in Te), 


It is clear that from a biological point of view 7 ...7°™ are 
much more realistic than 76? ...7°. Of the former T°!) seems 
to be the most realistic for reasons to be discussed elsewhere. 


Let us see now how 7‘® (P) maps onto the primordial graph P. 
The transformed graph contains m component primordials. In addi- 
tion each specialized point gives rise to v, subsidiary points. 
Each of the m primordial graphs receives 7, subsidiary points ac- 
cording to oe The total number of subsidiary points is, therefore, 


= i, . (16) 


Let n,, denote the total number of transport points in the primor- 


dial graph. In Figure 4, n,, = 9. Let n,,; denote the number of spe- 


cialized transport points in the ith component primordial. Then 


Me Shep i ees (17) 
i=1 


We now have 


n,=n-n,-(n,-7,,) . (18) 
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Because of (18), (17), and (12) expression (16) becomes 


tm 
>, F,=(n=n,) (m=1). (19) 
i=1 

Each of the (n~1,) (m- 1) subsidiary points is connected ac- 
cording to To to its own residual graph. Thus there are alto- 
gether m + as — 7.) (m—1) residual graphs. All their correspond- 
ing points and lines map on the corresponding points and lines of 
P in an m + (n—1,) (m — 1) to one mapping. 

All the subsidiary points belonging to a given specialized point 
map on the corresponding point of P. Every direct connection 
tf? —> f® maps upon the line f, — f, of P, while all the connec- 
tions (15) map on the zero connection "of a point to itself, in other 
words, on the point f, of P. The connections ine ——> ie 
(resp. f(?., <— f,)) which appear in 7° (2), as well as the con- 
TERS fe? we FP ER, esp. LO Fn KAO ed) 
and the corresponding connections of the eabsidiecy points map on 
the point f,, (resp. f/f. ;) of P. 

The expressions for the connectivity p of any of the transforma- 
tions 7°... 7? will, in general, be different from the one de- 
rived in I for T. We shall not discuss this now. 

Just as in I we prove the two theorems which hold for any of the 
fourteen transformations: 

Theorem 1. The transformation rules 7“ (u = 1,2,...,14) when 
applied to two different primordial graphs result, for the same 
choice of parameters, in different graphs. 

Theorem 2. The transformation rules T‘ (uw = 1,2,...,14) trans- 
form any partial graph of the primordial in such a way that points 
of 7‘ (P) which correspond to the points of P are connected in 
the same manner as those of P. 

Of greater practical importance is a theorem which justified the 
‘‘lumping’’ of several points of the primordial graph into one as 
discussed above (p. 118). 

Such a lumping may be considered as a special case of a many- 
to-one self-mapping of the primordial in which a partial graph P, 
of P is mapped on the point A of P; another partial graph P, of P 
is mapped on the point B of P, and all connections between any 
point of P, and any point of P, are mapped on the connection 
A—B. This can be done only if every line connecting a point of 
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P, with a point of P, is directed from the point of P,, to a point 
of P,. If there were some lines directed from a point of P, toa 
point of P,, then those lines could not be mapped on AB. 
Such a self-mapping, of which there may be very many different 
ones, results in a simplification of the actual primordial graph, re- 
ducing the number of its lines and points. A graph obtained by such 
a self-mapping will be called a simplified primordial graph P,. 
The graphs shown in Figure 1 of I and in Figure 4 of this paper are 
supposed to be simplified primordials. 

If we apply any transformation T‘?...7° to P, then the dif- 
ferent points of P,, unless they are all residual points, will in 
general be specialized in different component primordials of T(P). 
The same is true about the points of P,. Each of those special- 
ized points will, in general, give rise to several subsidiary points. 
Because of oe (u=1,...,14) and FY some of those points, 
which we shall call Q-points and which arise as a result of spe- 
cialization of points of P,, will be connected to some points, 
R-points, which arise as a result of specialization of points of P,. 
However, because of (es and TS, every point Q of T™ (P) is 
connected directly or indirectly with a point R of T™ (P) in such 
a manner that the arrow leads from Q to R: Q —> R. 

Denote by 7“” (a) the set of points of 7‘ (P) which arises as 
a result of specialization of the point a of P under the transforma- 
tion T™, Then 7“ (A) is a subset of the set of the Q-points, 
and 7‘) (B) is a subset of the set of the R-points. We can now 
map in several different ways the set Q on T“ (A) and the set R 
on 7“ (B), mapping each point of T“ (A) and of T™ (B) onto it- 
self. The residual graphs belonging to the Q-points and to the 
F-points are correspondingly mapped on the residual graphs of the 
point of 7™ (A) and T™ (B). As a result of such a mapping we 
obtain a simplified transformed graph which is the transformed 
graph of the simplified primordial P,. Thus we have 

Theorem 3. There exists a (rs (P)], which is identical with 
ee (P ey: 

All empirically determined physiological graphs are simplified 
graphs. Theorem 3 enables us to operate with simplified graphs 
and to obtain conclusions which can be verified experimentally. 
Thus in a simplified primordial graph absorption f, is connected 
via f,;, to synthesis f,, which is a simplification. In a trans- 
formed graph, for a given choice of parameters, there will be con- 
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siderably more complex connections between groups of absorption 
phenomena and groups of synthetic processes, each of which rep- 
resents a ‘‘lumping’’ of a still greater number of such processes. 
The simplified transformed graph may be still amenable to experi- 
mental study, whereas the actual complete transformed graph may 
be too complicated for it. It may be possible by considering pri- 
mordial graphs of a different degree of simplification to work out a 
method of successive approximations for the study of the trans- 
formed graphs. 

It seems likely that theorem 2 of III with the correction given 
subsequently by E. Trucco (1956), as well as the theorems of IV, 
hold for all transformations JT‘? ...7°0”. We shall only outline 
here the reasoning which would be used in studying this question. 

For transformations T‘® ... 7 we cannot consider that all 
corresponding specialized points are merely fused. This holds 
only for such specialized points which are not adjacent to a trans- 
port point alte (1); (w=8...14)]. For other points this does not 
hold because of 7§” (2) and T{”(3). In fact, while in the trans- 
formations T°” ...7™ a specialized point did topologically be- 
long to all m component primordial graphs, though biologically it 
belongs to only one, in 7“ ...7“ each specialized point does 
topologically belong to only one given component primordial. This 
is another advantage of the transformation T7‘®...7°. Never- 
theless any point of a component primordial graph z in T™ (P) (u= 
8...14) can be reached by a directed way from Be , where BS” is 
defined on page 121 of III. 

The ith component primordial graph is obtained from P by remov- 
ing n, points. Consider a point p™ in the ith component primor- 
dial. If the corresponding point p of P can be reached in P from 
any base point by a directed way which does not pass through one 
of the 7, removed points, then our assertion is clear. Let, how- 
ever, p be reached in P from a base point only through a point g 
which is one of the 7, removed points. Let the directed path from 
a base point a to p be 


a—>b—->c—>q->d—gj—p. (20) 
The points 


a®, 6, eo, d®, g®, and p@ 


are present in the zth component primordial, but gis not. Either c 
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or d or both must be transport points. In that case we can reach 
the point p‘ from the point a‘ by the following directed path, 
which exists because of (nD (2) and EO3) (U8 oi 4): 

When the transport points are not specialized: 


: ; : k k) (k) (¢) 
a) —, BO 5 f% > (erg ate EE ree ae 


; 21 

g® oats p@ ; (21) 

When the transport is specialized in the sth component primordial: 
: ; , k k 

a ao be ae peer ee oe Pas ee sae q‘ } reed Pee as 


) OMe 


fetes alr heme 
Similar considerations can be used in the proofs of the theorems of 
IV. 

All this suggests the possibility of formulating general condi- 
tions which a transformation should satisfy, without specifying the 
transformation, in order that all the theorems of III and IV should 
hold. 

Equation (1) of IV is now somewhat modified due to the fact that 
in T® we consider nm, instead of n — n, and in general the inequality 
sign in (13) holds. 


SY 


The author is indebted to Dr. Ernesto Trucco for a discussion of 
this paper. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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The role played by the group of a graph in determining sets of equiva- 
lent points is exposed and illustrated by a few simple examples. 


N. Rashevsky (1955) has recently introduced the concept of top- 
ological information content of a graph in view of possible biologi- 
cal applications. 

Any graph may, in fact, be used to realize a discrete source of 
information according to the abstract definition given by Shannon 
and Weaver (1949, p. 9). This is done by considering each set (or 
‘*class’’) of topologically equivalent points in the graph as a ‘‘let- 

’ of an ‘‘alphabet’’ and assigning to it the probability 


OF : 
eh ial (1) 
where the index j enumerates the sets and n, is the number of 
points in the jth set. We have 


y= (2) 


m being the total number of sets. The quantity w, represents the 
probability of selecting a point belonging to the jth set in a random 
choice and may be called the weight of that set. In sending mes- 
sages, each ‘“‘letter’’ is chosen with probability w,. If these prob- 
abilities are independent, i.e., if the choice of ae one letter is 
not influenced by the previous choices, we define the information 


129 
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per symbol (also called information content or entropy of the 
source) by the equation 


[=- Ne wlgw, bb (3) 
=i 


It is the purpose of this note merely to point out the relation be- 
tween the idea of topological equivalence on one hand and the 
group of the graph, as defined by D. Kénig (1936, p. 5), on the 
other hand. It should be noticed that this, as well as several other 
pertinent problems, is discussed in a very long and extremely in- 
teresting paper by G. Pélya (1937). 

All the one-to-one mappings of a graph onto itself form a group 
which is called the group of inner automorphisms or simply the 
group of the graph. It is clear from the work of N. Rashevsky (doc. 
cit.) that a point of the graph can only be transformed into a topo- 
logically equivalent point by an operation of the group. Accord- 
ingly, we may use the following 

Definition: Let H be the group of the graph. Two points, say 
a and 8, of the graph are equivalent if and only if there exists at 
least one operation of H which maps a onto 6. Denoting this equiv- 
alence by a~b, we see that the above definition satisfies the 
three basic requirements of equivalence, namely, 


1 acer: 
2) liao bob <a = 
3) Ifa~ b, b~c, thena~c. 


The situation is similar to that discussed by Speiser in connec- 

tion with the symmetry of ornaments (cf. Speiser, 1937, $31, 1st 

ec., p. 95), There, however, the relevant group is that formed by 
the translations and rotations of plane figures. 

Since the elements of a graph are points and lines, its group 
can be written as a permutation group of points and lines, and 
incidentally, the concept of equivalent lines can be introduced in 
a similar way as was done for the points. Consider, as an ex- 
ample, the very simple graph of Figure 1. Its group H consists of 
the four operations F (the unit), A, B, and C, which can be speci- 
fied as follows: 


*We use the following notation: 1g—for logarithms to the base 2, 


log—for ordinary logarithms (base 10), and In—for natural logarithins 
(base e), — 
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E: (a) (0) (c) (d) () (B) (y) (8) (6) , 
A: (a) (0) (ed) (x) (By) (8) (©) , (4) 
B: (a) (8) (c) (d) (&) (B) (y) (8e) , 
C: (a) (0) (ed) (a) (By) (8e) , 
and satisfy the group relations 
AA=B*=C7=£, 
AB=BA=C. 

In (4) the permutations are written in cyclic notation; e.g., (cd) 
(By) means that the point c is to be mapped onto the position of 
d, d onto c, and similarly for the lines 8, y. The group (4) is 
clearly the direct product of two groups of order 2. Normally, 
cycles containing only one element would not be written down. 
Also, the permutations of the lines need not be written if the graph 


is such that not more than one line joins any two points. (Frucht, 
1938, $1, p. 240). 


ed 
a 
eb 
B + 
r) 
Ce ed 
ee er 
ic 
FIGURE 1 


It was pointed out by Pélya (loc. cit., p. 158) and follows im- 
mediately from our definition of equivalence that the sets of equiv- 
alent points (or lines) in the graph are nothing else but the ¢ransz- 
tive sets of elements in the group H (e.g., cf. Speiser, loc. cit., 
§36, p. 110; or Burnside, 1911, $133, p. 170; $142, p. 185).* 

A theorem due to Frobenius (cf. Burnside, loc. cit., pp. 189-190) 
may be useful in this connection. If n is the degree of the group 


*For this reason we prefer to speak of sets of equivalent points, and 
not classes, as was done by N. Rashevsky (loc. cit.) and later by G. 
Karreman (1955). The word ‘‘class’’ is usually reserved for classes of 
conjugate operations in the group. 
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(i.e., the number of elements permuted) and & its order, and if, 
furthermore, ¢ is the number of transitive sets and v, the number of 
operations of the group which leave exactly r elements unchanged 
(r= 0,1, 2,..., 7), then we have 


n 


ye vi=h 


r=0 


(this is trivial), and 


ye rvi=ht. (5) 
rol 


Since two conjugate permutations are similar, i.e., have the 
same cycle structure (Burnside, loc. cit., $10, p. 7), it follows 
that the sum in equation (5) can have at most p non-zero terms, 
where p is the number of classes of conjugate operations in the 
group. If the operations of the group are represented as permuta- 
tion matrices of n rows and columns, the numbers 7 are the char- 
acters of the matrices in this (reducible) representation. (Speiser, 
loc. cit., $389, p. 118; also cf. Pélya, loc. cit., sec. 20, p. 167). 


By using formal power series Pélya gives a general theorem 
enabling one to find the number of non-equivalent configurations 
if the group H is known. 

The definition of equivalence given above remains valid if H is 
any group which permutes the points of the graph. For certain 
problems it may be useful to define the group H as the totality of 
mappings which transform a graph into a congruent graph. This 
implies that one should give a careful definition of congruence be- 
tween graphs. This becomes particularly important when we are 
dealing with graphs in which the points may be physically dis- 
tinguishable (e.g., each point the seat of an atom) or, as we shall 
call them, differentiated graphs. The congruence may be defined 
in several ways (Pélya, loc. cit., secs. 32-35, pp. 183 ff); each of 
them has a definite significance for chemical applications. De- 
pending on which definition of congruence we choose, we obtain a 
different group of the graph; this, too, is illustrated very clearly 
by an example in Pélya’s paper (sec. 56, p. 212). 

Similar considerations will apply to the lines if they become 
differentiated, e.g., oriented. 
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These remarks should be useful when one attempts to calculate 
the information content of some actual simple molecules, or rather 
a family of isomeric molecules, as was recently done by G. Karre- 
man (unpublished), 

In the following we give two simple examples showing the rela- 
tion between the group of the graph and its sets of equivalent 
points. We shall consider graphs to be congruent in what Pélya 
calls the ‘‘topological’’ sense, i.e., we use his conditions (1), (Il), 
and (III), but not condition (IV) (Jd2zd., p. 186). 


Ce eda 


Senora 


o Ee 4 b 
= Wy; *B | | 
c d 
FIGURE 3 
eS 
FIGURE 2 


a) For the rectangle shown in Figure 2, it is sufficient to notice 
that its group of symmetry is transitive. It consists of the four 
operations 


E 
A = (ab) (ed) , 
B= (ac) (dd) ’ 


C = AB = (ad) (bc). 


The group of the graph must therefore be transitive, and there is 
only one set of topologically equivalent points. 
b) The graph of Figure 3 has the following group of order 8: 


E 
A = (ad) ’ 
B =(cd) , 


C = (ac) (bd) (9) , 

D = AB = BA = (ab) (cd) , 
F = AC = CB =((adbe) (fg) , 
G = BC = CA = (acbd) (fg) , 
H = ABC = (ad) (bc) (f9) . 
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Equation (5) is readily verified, since h=8, n=9, ¢= 5, 
5 1 =¥y =U, =e =, = 9,43 =4, 0, = 1, ¥, = 2, vy = 1. The 
transitive sets of points are 


[a, 2, ¢, 2), {fg \, te), (4, Let. 


All the examples given in N. Rashevsky’s Figure 1 (loc. czt.) 
can be treated in a similar way. For instance: 

For Figures 1b, 1d, 1g, 1i, and 1k, the group H consists simply 
of the unit element E. (This implies that congruent graphs must 
have their lines oriented in the same manner.) Hence, / =lgn. 

For Figure 11, the group H consists of two elements: FE and A = 
(1) (25) (84) (67). Therefore 


Vv, =v 


I : lg 7 : l d 
ne See ter 

For Figures la, le, and 1f, the group H is transitive, and thus 
T=. 

In the same way we can determine the groups of the graphs given 
in Karreman’s paper (loc. cit.). 

Moreover, Rashevsky’s theorems 1 to 5 can be proved by finding 
the groups of the corresponding graphs. 

We may also mention the following rather trivial theorem: Let G 
be a graph; to each point p of G assign a supplementary point p’; 
join p and p’ by a non-directed line. Denote the graph thus ob- 
tained by G’. Then, neglecting the contributions due to the lines, 
the topological information contents of G and G’ differ by one bit. 
Except when G consists of one point only, the graphs G and G’ 
have the same abstract group (Frucht, loc. cit., $2, p. 241), but 
clearly to each set of topologically equivalent points in G there 
correspond two sets of G’ with the same number of points. Hence, 
if n and m refer to G: 


These considerations are not likely to be very useful in practice 
Since to find the group of a graph (except for the simplest ones) 
is a rather difficult problem. (For a special case see: Frucht, 
1937; for non-directed trees: Pélya, loc. cit., section 54, p. 208). 


Some interesting remarks on the groups of molecules, their rela- 
tion to isomerism, and the use of Frobenius’ theorem to obtain the 
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number of transitive sets can be found in Lunn and Senior’s paper 
(quoted by Pdlya). 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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MATHEMATICAL BIOPHYSICS OF COLOR VISION 
Il: THEORY OF COLOR CHANGES INDUCED BY THE 
ALTERNATION OF COLORS AT VARIOUS FREQUENCIES. 
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A model previously introduced to account for a number of phenomena 
of color vision is studied for the conditions in which the alternation of 
colors is periodic. The case of three primary receptors having no dif- 
ferences in their time constants is considered for the situations such as 
the temporal alternation of two primaries, a primary and a neighboring 
binary, a primary and white, a binary and white, etc. Using the same 
mathematical model which was used to account for the enhancement ef- 
fect of interrupted illumination, it is found possible to account for ob- 
served changes in hue in a qualitative manner. A method is suggested 
for measuring the color changes quantitatively. This method is readily 
adapted for demonstration purposes. The model shows that there should 
be differences between ‘‘primaries’’ and ‘‘binaries’’ with respect to 
change in hue. In principle, therefore, comparison between theory and 
experiment should yield information regarding the physiological 
primaries. 


In a previous paper (Landahl, 1952), familiarity with which is 
presupposed, a mechanism was proposed by which the excitation 
of a small number of primary receptors each having its own sensi- 
tivity curve for various wave lengths of light results in the ex- 
citation in a certain region of a higher center C’, the position of 
this region depending upon the extent to which the primaries have 
been excited. A number of the phenomena of color mixture can be 
accounted for by this model. These are the phenomena which arise 
in successive contrast. We shall consider only steady-state ef- 
fects, and hence the experimental situation is that in which colors 


*This research was supported by the United States Air Force through 
‘the Air Force Office of Scientific Research of the Air Research and 
Development Command under contract No. AF 18 (600) 1454. 
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are presented in a periodic manner at constant frequency, such as 
with a Maxwell disc. In the steady state only the dominant time 
factors are apt to be important so that the model we introduce does 
not need to be so complex as to account for all the details of after- 
images. It is the purpose of the present paper to consider some of 
the effects which result when the color mixing occurs at frequen- 
cies which are so rapid that the excitatory processes cannot follow 
the changes accurately but not rapid enough for the excitatory 
processes to respond to mean values only. 

Thus we may hope to account for some of the effects which oc- 
cur when colors are presented in various temporal patterns. It may 
also be possible to estimate some of the parameters of the system 
and obtain information regarding the physiological primaries. 

For the present we restrict ourselves to a very special case, 
only attempting to illustrate some of the effects which can be ac- 
counted for in terms of the model. In order to simplify matters we 
shall suppose that there are three primary receptors. The param- 
eters of all the receptors are the same except for those param- 
eters which determine sensitivity to wave length. Furthermore, 
we shall use the model in which the excitatory and inhibitory proc- 
esses are linear in the stimulus intensity S (/bzd.). Let a stimulus, 
having a given spectral energy distribution, be characterized by 
the quantities S,, S,, and S, which represent the effective stimuli 
acting upon the three respective receptors, the units of S being 
chosen in such a manner that S, = 8, = 8, for white light. We 
Shall suppose that when the stimulus S, is initiated there develops 
an excitatory process ¢, and an inhibitory process j, at the ends of 
the neuroelements corresponding to the ¢th primary. The quantities 
€, and j, are assumed to be governed by the equations (Householder 
and Landahl, 1945, chap. i) 


de. 
a = a(¢, > eo ’ (1) 
dj. 
= ee b(R¢, Sa i;) ’ (2) 


in which the rate constants a and 6 as well as the ratio R (= ¢/y, 
Tbid., chap. vi) do not have subscripts since they have been as- 
sumed to be the same for all z. The quantities @, are assumed to 


be linear functions of S,, the constants being assumed the same 
for each value of z. 
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We shall restrict our considerations to the steady-state situ- 
ation when two colors are presented alternately at a given fre- 
quency. Let 7 be the frequency, equal to 1/7, T being the pe- 
riod. If during a fraction r of the period the intensity is such 
that ¢,=¢%,,;, but during the remaining time (1-1) T the inten- 
sity is such that &;=%,;, we can determine ¢, from (1) as a 
function of time in the interval from 0 to r7’, using ¢ = o,- Then 
using ¢ = $ a» together with the value of ¢,(77) already found, 
we can find «, as a function of time from rT to T. Continuing 
this process for the intervals T to 27, etc., we find a solution 
which becomes more and more independent of the initial condi- 
tions and depends only on the instant ¢ in the cycle. The de- 
sired solution in the steady state can be obtained as a slight 
generalization of results by N. Rashevsky (1948 chap. xxxvii; also 
cf. Householder and Landahl, chap. vi). If we introduce 7= (1-7), 
the result is as follows, for ¢ = 1, 2, 3: 


O< tS rT, 
: 7 e~ ar i ee e~ ar 
—at at 
sf Pai ! 1 aT © a d 2; Bry ? 
e l-—e 
(3) 
mr StST, 
ees e7 ur es eur 
~a(t—rT) ~a(t-rT ) 
c= bau oar P| py [1 ar ‘ 


The expressions for j, are obtained from (3) by replacing a by 6 
and ¢ by Rd. The significant quantity is the difference between 
the excitatory and inhibitory factors which we shall denote by a, , 


These quantities, 7, , 7,, and o,, will determine wnich region in 
the center C’ (Landahl, 1952) will become excited. This will 
always fall within the triangle C7C3C;-. If Q° is the center of 
the active region, corresponding to a stimulus, having the com- 
ponents o, , 9,, and a,, then its position can be characterized by 
the quantities OX’, OY’, and OZ’ of Figure 1. It can be shown 
that OX’ + OY’ + OZ’-=0, and only two values are needed to de- 
termine the position. We shall hereafter drop all primes. Then, if 
X is equal to the distance OX, Y = OY, etc., we have the position 
of the representative point, given as follows in terms of 0,, %, 
and a, [/béd., equation (6) with y; = ¥2 = Ys }: 
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é 20, -0,- 93 (5) 

2(0, + % + 95) 

when all o’s are positive. If equation (4) yields negative values 
for any o, , then in (5) we put o, equal to 0. Note that X is a func- 
tion of position ¢ in the cycle through the o’s. In (5) we have also 
ignored the thresholds and hence we must consider situations in 
which the intensity is relatively high. Corresponding to (5) are 
expressions for the components Y and Z along OC, and OC, , ob- 
tained by permuting the subscripts, i.e., changing 1 to 2, 2 to 3, 
and 3 to 1. If the expressions for Y and Z ate written out, it be- 
comes evident that X + Y + Z = 0 as mentioned above. 


FIGURE 1 


Before continuing with the quantitative model we shall illustrate 
graphically how changes in color may be brought about in this 
model. Consider a situation in which a primary color is alternated 
with white, so that S, is a constant, s? , and S, = S, = 0 in the 
interval from 0 to st; while §, = S, = S, = 35 from ae to T. From the 
nature of the equations, which can account for the Briiche-Bartley 
phenomena (Bartley, 1951), the net excitations 0,,%,ando, will 
vary in time as illustrated in Figure 2. The stimulus pattern is 
shown in the same figure. It will be noted that, during most of the 
time from ¢ = 0 to ¢ = T/2, the primary color will appear somewhat 
more intense. (If an unsaturated primary were used, it would also 
appear more saturated.) From —T to T, following a large reduction 
in the intensity of S,, there is a decrease in o, resulting in its 
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average value being less than its steady-state value. On the other 
hand o, and o, tend to overshoot during the second half of the 
cycle. As a result during this time the values of ,,%, ando, 
are no longer equal, but o, =o, >o, for most of the interval. 
Hence the white fraction appears as a color which is complemen- 
tary to the primary color. 


0 Yat T oh geen ena T 
FIGURE 2 


In order to illustrate what phenomena can be accounted for by 
this model, we shall give the results of some numerical calcula- 
tions. For this we choose a/b =10 and R =— and divide each 
period into 10 equal intervals, calculating the position of the re- 
gion stimulated at the ends of those time intervals. The position 
of this active region in the center (',C,C, determines the color 
experienced. The numerical values of the parameters are chosen 
for convenience of calculation, though they show the qualitative 
features of the response of the ‘‘continuous’’ fibers of the retina 
as well as of the enhancement effect in the alternation of black 
and white (Bartley, loc. cit.). 

Case 1. Alternation of two primaries. Consider the case in 
which a primary color P,, corresponding to the point C, where 
X=1, Y=Z=~-4, is alternately presented with the color P,, 
corresponding to C, where Y=1. Let the intensities be equal. 
_ Then for the time ¢ in the interval 0 < ¢ < - T, $4, =%y> and d,, = 
$a, = 0; while for +T<t<T, $,, =, =9, and dm =o. The 
results for f= 0, 1, 5, 10, 20, 40, and ~ (in units of 6) are shown 
in Figure 3. Since ¢, is always zero so that Z = - 3, the 
region stimulated remains near the line C,C,. Hence only this 
line is shown for each of the frequencies, the length being equal 
to /3 if the distance OC, is used as the unit. When the frequency 
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approaches zero, the point representing the active region simply 
jumps back and forth between the points C, and C, corresponding 
to the primary colors used, only a negligible amount of time being 
spent between these points. When f= 1 and ¢ = 0.17, the point is 
a distance 0.2 from C, while from 0.27 to 0.5T the point is at C,. 
By 0.67 it has passed all the way to within 0.2 of C,, remaining 
at C, for the rest of the cycle. Hence there is essentially an al- 
ternation between the primaries P, and P,. At a frequency 5, the 
representative point starts at the time 0.17 at 1.22 (number 1, 
Figure 3), a point somewhat closer to the mid-point than to C,, 
and moves rapidly toward C,, almost approaching it at ¢=0.5T 
and then slows down somewhat. But after this it begins to move 
rapidly back toward C,. In this case the points in Figure 3 are 
numbered in order, return values being indicated by marks below 
the base line. At higher frequencies the range decreases until by 
f = 40 only a small region about the mid-point is active. 


C, # 4 6 .8(C;G) 1.0 EZ 1.4 16 Co 


fe gee Sees 3 2 oh 
6 7 8 9 10 


FIGURE 3 


If, for example, we assume that red and green are primaries then 
at very slow speeds there should appear a simple alternation of 
colors. At a somewhat higher frequency the red and green should 
tend toward orange and yellow-green respectively; the possibility of 
the appearance of yellow would be low. At a still higher frequency 
the red and green should become orange and yellow green, with 
yellow possibly appearing between. Note that because of the rapid 
motion of the representative point, the maximum excitation tends 
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to decrease and only at higher frequencies does it increase as the 
representative point approaches its final position. It seems plau- 
sible to suppose that the certainty of the color recognition is re- 
lated to this intensity. Then at these intermediate frequencies 
(e.g., f= 10 or 20) the recognition of the colors will be poor as 
compared to that at higher frequencies. 

A quantitative estimate of such an intensity can be fairly easily 
made for the mid-point and for the ends if we suppose that, at any 
point on the next higher center, C, is represented by neuroele- 
ments of the simple excitatory type with a rate constant, a, the 
same as that of the primaries. Then if we approximate the time 
course of excitation in the region by a square wave of diameter 
A on the triangle, or width A along C,C, (Landahl, 1952), then we 
may use expression (3) to estimate the excitation at any time. 
Here $3; = 0, and ¢,,; = $¢’, a constant which we do not need to 
know. The ¢’s represent the mean values of the net excitation 
in the region A. For the mid-point, 7 must be replaced by 7/2 and 
rT by the time it takes the wave of width A to cross the mid-point. 
This can be estimated from the Figure 3. For the estimation of 
the intensity at a distance A from the end point, T is not changed 
and r7 is replaced by the total time any part of the square wave 
remains at this point. This also can be estimated from data given 
in the figure together with some additional points. If A = 0.1, then 
we find at the mid-point between C, and (, the following relative 
excitation values for the various frequencies: (0.07, 1), (0.07, 5), 
(0.17, 10), (0.24, 20), (0.45, 40), (1, ~); the relative values are 
given first and the frequency second. On the other hand, near the 
point of maximum excursion toward C, or C, we find the following 
values: (0.99, 1), (0.18, 5), (0.15, 10), (0.29, 20), (0.41, 40). 

Case 2. Alternation of a primary with white light of equivalent 
intensity. Consider next the case in vnach a primary, Heiss is pre- 
sented with intensity / from the time 0 to > 17, while from + : —T toT 
all primaries are equally stimulated with an intensity 4+/. Using 
the same values for the parameters as in case 1, we obtain the 
following results. At very slow frequencies ihe representative 
point remains at either C, corresponding to P, or at the origin, 
corresponding to white (W), except for a Reais time spent in 
transit. At somewhat greater frequencies (e.g., f= 1 in units of d, 
Figure 4) there is an alternation between C, and a point corre- 
sponding to an unsaturated color, the complement of P,. At higher 
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frequencies the recognition of color may become uncertain; P, and 
its complement are absent. The point corresponding to the ex- 
cited region moves back and forth along the line OP, but it no 
longer reaches either P, or O. However, the likelihood of re- 
sponse, as discussed under case 1, is somewhat greater near the 
points of maximum excursion along OP,. But this is no longer the 
case if the frequency is increased too much. At considerably 
greater frequencies, the region of response narrows down toward 
its final value, the certainty of that response increasing with 
frequency. 


f =00 SS ee eee 
q W (Co C3) 
FIGURE 4 


If instead of giving equal times to P, and white the proportions 
are — and 2, then we find for the values of the parameters used 
that, while the extent of the saturation of the complement of P, is 
less than in the previous case CG xy EL the probabilities of the re- 
sponse to it are somewhat increased. This may be related to the 
fact that under some conditions the primary is not seen although 
its complement is seen at a low saturation. 

Case 3. Alternation of a binary with white. The result of al- 
ternating a binary with white is illustrated in Figure 5. It can be 
seen that in this case the saturation of the complement is con- 
siderably less than that shown in Figure 4. This suggests a 
method for obtaining information on the physiological peleiee. 

Case 4. Alternation of binaries. Let a binary, P ,P, (ie., — =P, + 
—Ps)s be alternated with another binary, P,P,. At. vom low apa 
the representative point jumps back and ‘torth between the point 
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C,C, (midway between C, and C,) and the point C,C, (Figure 6). 
When the frequency f is equal to 1 in units of 6, we find that the 
alternation is now between a first point near P,P, but in the di- 
rection of P,, and a second point between C,C, and C, but much 
nearer to C,C,. This effect can only occur in this model if P,P, 
and P,P, are both fairly close to maximum saturation. 

Case 5, Alternation of a primary with a neighboring binary. The 
results of alternating a primary (P,) with a neighboring binary 
(P,P,) are illustrated in Figure 7. Here it is seen that at a low 
frequency (f= 1) there is essentially an alternation between the 
primary (P,) and the color lying between the original binary (P ,P,.) 
and the primary P,. Thus the apparent binary moves away from the 
primary. At higher frequencies this effect disappears and the 
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FIGURE 6 


146 H. D. LANDAHL 


representative point becomes spread out over a region which, with 
higher frequencies, condenses toward the point midway between 
C, and C,C,. If we suppose that yellow is a binary (P,P,) or at 
least has a component which is a binary (Hartridge, 1950, p. 279) 
and that P, is red and P, is green, then the above can account for 
the appearance of a greenish yellow when red is alternated with 
yellow at a low frequency. Similarly yellow takes on an orange 
tint when alternated with green, and a yellow orange may appear 
more yellow if alternated with red. 


C, 2 yi & 8 1st 125 gh et eee 

ooo ooee—ae—oeee—oe TTT 
(Z=1/2) 

SCALE (C,C) V; 
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FIGURE 7 


Case 6. Alternation of white or colors with black. Alternation 
of white and black leads only to transient changes in intensity 
discussed previously (Householder and Landahl, 1945, p. 45). This 
occurs in the model because of the assumption of equivalence 
made for simplification. If this assumption is dropped then, in 
principle, transient colors could appear if the time constants (a,) 
are sufficiently dissimilar (cf. Piéron, 1943). Because of the ef- 
fect discussed under case 1, the differences must be quite marked 
in order that colors would appear even under optimum conditions. 
It may be noted that the Fechner Type I induced colors (cf. 
Linksz, 1952) are markedly influenced by intensity and tint of the 
illumination, the colors being only appreciable when the illumina- 
tion is low. On the contrary, the steady state successive con- 
trast effects mentioned above are relatively independent of in- 
tensity. In the case of Fechner Type II colors (Benham disc) 
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spatial factors enter (simultaneous contrast). We shall not con- 
sider this problem here. 

When any color is alternated with black at low frequencies the 
only result under the assumptions of the present model is that of 
intensity changes. Even though such transient changes are diffi- 
cult to measure, from the point of view of the present theory such 
measurements would be useful because they would help to de- 
termine the parameters of the system. 

The color changes due to successive contrast in the steady 
State can be most readily demonstrated with discs in which colors 
A and B are sectors of desired angular widths. This has been 
done in some preliminary experiments which confirmed the above 
theoretical conclusions. For comparison, a central area or peri- 
pheral ring of color A or B can be used. For the influence of color 
A on B, the sector A may be made two or three times the angular 
width of B, the central area being of color B for comparison. The 
two areas should be separated. A gray ring separating the central 
area from the rest of the disc tends to prevent the sector of color 
B from appearing to be part of the central disc. Having demon- 
strated a change in the color of B, the central area may then be 
changed to match the apparent color of the sector B. The central 
area may be a Maxwell disc rotating at a rate necessary to produce 
a uniform color which may serve as a standard for comparison. 
Thus for example, red and yellow, in the ratio two to one, rotated 
at a rate of a few revolutions per second give the appearance of 
red being alternated with a greenish yellow. The latter can be 
matched by rapidly rotating two parts of yellow to one part of 
green. The shift of yellow toward both green and red may be 
demonstrated simultaneously by adding to the above disc an outer 
ring having green and yellow with angular widths in the proportion 
two to one. A narrow gray ring should be used to separate the new 
ring from the one just central to it. When such a disc is rotated 
at a rate of a few cycles per second the contrast between the ap- 
pearances of the two yellow sectors becomes very marked. A 
similar situation obtains if yellow is replaced by purple and the 
green by blue; or if yellow is replaced by a blue-green and the red 
by blue. 

It has been pointed out above that primaries which are quite 
saturated should be very resistant to change in hue. Some pre- 
liminary qualitative observations show that reds, greens, and blues 
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show some increase in saturation and an increase in intensity as 
expected, but little change in hue. 

If one of the colors A or B is white (or gray) and if that sector 
is fairly small then there is a change in hue which depends upon 
the illumination. Under low illumination and under yellowish light 
(ordinary lighting), a bright sector on a dark background appears to 
have a yellow or pink hue. This is very pronounced with white on 
black. This effect can be avoided by using moderately strong il- 
lumination of daylight whiteness and by using sectors which are 
more than forty five degrees. 

Since there is considerable evidence of more that three pri- 
maries, (Hartridge, loc. cit.; Granit, 1955), it would seem desirable 
to consider the situation in which some or all of the ‘‘binaries”’ 
have some primary properties. The interesting results of K. Moto- 
kawa and M. Ebe (1952) should also be taken into account. On the 
basis of the present model, the presence of additional types of 
receptors does not necessarily alter the basic theory requiring 
three colors for matching a given color, except in regard to satura- 
tion. However, the changes in hue, saturation, and intensity being 
discussed are dependent upon the properties of the actual recep- 
tors. Thus, for example, suppose that yellow is both a primary 
receptor and a binary in the sense that it can be stimulated by red 
and green independently of yellow. Let p,(A) be the sensitivity 
as a function of wave length for the yellow receptors which have 
a maximum sensitivity in the yellow. Let p}(A) be the sensitivity 
of the receptors having one maximum in the red and one in the 
green. Then two spectral distributions of light which give the 
same apparent color may be different in their changes when al- 


ternated with other colors or may have different effects upon other 
colors. 
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Experimental data concerning stress-strain relations of biological sub- 
stances, including single muscle fibers, are reviewed in the literature. A 
reconsideration of the conclusions reached from these data is made. 
Thus, the experiments indicate that prima facie a linear theory of elas- 
ticity such as Hooke’s law does not obtain because of the complicated 
form of the stress-strain curves. It is then shown that from a considera- 
tion of true stress at a point, where true stress is the force per unit ac- 
tual area at the point, Hooke’s law may be so generalized as to predict 
qualitatively and semi-quantitatively the complicated experimental curves. 
This complication arises from the long range ‘‘elasticity’’ of biological 
substances in general and muscle fibers in particular. Thus, substituting 
true stress for stress, the stress-strain relation of muscle fibers remains 
linear. In the case of ‘‘small’’ strains, the true stress reduces to the 
usual stress concept in which the force is considered per unit initial area 
at the point. 


I. General. 


In recent years, a great deal of interest has centered about the 
mechanics of muscle fibers and, in particular, the complex tension- 
length diagrams to which they give rise. It is apparently an a 
posteriori conclusion among most of the investigators that except- 
ing small stretches Hooke’s law is, in general, not obeyed (Ramsey 
& Street, 1940; Guth, 1946; King and Lawton, 1950; Buchthal, 
1942). 

Experiments on myosin filaments have similarly given rise to 
nonlinear curves (Dubuisson and Monnier, 1943). These curves are 
often difficult to interpret on physiologic grounds and give rise to 
many perplexities (Ramsey, 1946). In an attempt to account for the 
nonlinear curves produced by these biological substances, Bull 
(1946) applied elastomeric theory to Ramsey and Street’s data on 
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single muscle fibers, but could not account for the shape of the 
curve in the range of large extensions. In the present communica- 
tion, an attempt is made to obtain curves similar to that of Ramsey 
and Street by the use of Hooke’s law. Although the derivation of a 
nonlinear stress-strain relationship from a linear one seems at first 
glance to be contradictory, it will be seen that the apparent para- 
dox rests on a methodological procedure involving the ‘‘form’’ of 
Hooke’s law in the case of large deformations. 


NORMAL FORCE-LENGTH DIAGRAM OF SINGLE FIBERS 
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FIGURE 1. Length-tension curves for active (stimulated) muscle fiber. 
Ordinate: force per unit cross-sectional area in arbitrary units. Abscissa: 
length of muscle in per cent of resting length. Solid line calculated to 
fit data, from equation 3, of R. W. Ramsey, Ann. N. Y. Acad. Science; 
1946, 47, 675. (Reproduced by permission from original of figure 4, R. 
W. Ramsey, Ann. N. Y. Acad. Science, 47, 681). 


In carrying out the analysis which appears below, reference is 
made to the fact that the definition of elasticity is independent of 
the magnitude of deformation at a local point of the medium (Planck, 
1932), ideal elasticity being concerned only with the return to equi- 
librium of the displaced points after the load is removed. The mag- 
nitude of the strain or elongation is not restricted, i.e., the stress 
is everywhere proportional to the local deformation (Hooke’s law). 
The actual elasticity of all substances is, of course, modified by 
the previous history of treatment, hysteresis, and so forth. All sub- 
stances are elastic within a certain range. In general, metals have 
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exceedingly small ranges, whereas rubber retains its elastic prop- 
erties under relatively large deformations (Novozhilov, 1953). Bio- 


logical substances exhibit the latter property (King and Lawton, 
loc. cit.). 


Il, The Form of the Tension-Length Curve. 


Reference to the literature (Ramsey & Street, oc. cit.; Buchthal, 
loc. cit.; Dubuisson and Monnier, loc. cit.) in most cases shows 
the form of the tension-length curve to be nonlinear. Figures 1 and 
2 from Ramsey & Street reproduce the curves for stimulated and rest- 


8 


100 120 140 160 180 


FIGURE 2. Length-resting tension curves. Abscissa: length in per cent 
of resting length. Ordinate: tension in arbitrary units. (Reproduced by 
permission from original of figure 2, R. W. Ramsey, S. F. Street, J. Cell 
and Comp. Physiol., 15, 18). 


ing muscle, whereas Figure 3 is the total tension-length diagram, 
over the entire range. Guth’s analysis (oc. cit.), based on elas- 
tomere theory, attempts to explain Ramsey & Street’s data by com- 
parison with that for rubber (Figure 4). The range of strains 
0 <e«<e* corresponds to the tension diagram of active muscle 
(Figure 1), whereas that for ¢ > e* corresponds to the curve of rest- 
ing tension (Figure 2). The total tension-length curve (Figure 3) 
is then approximated by the rubber stress-strain curve (right-hand 
curve of Figure 4). 


Ill. Kinematical Consideration. 


If dr is the vector distance between two particles before the de- 
formation (translation and rigid body rotation excluded) and if the 
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FIGURE 3. Length-total tension diagram of an isolated muscle fiber. 
Ordinate: total tension in per cent of maximum at resting length. Ab- 
scissa: lengths in per cent of the resting length (100). (Reproduced by 
permission from original of figure 4, R. W. Ramsey, S. F. Street, J. Cell. 
and Comp. Physiol., 15, 21). 
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FIGURE: 4. Left: stress-strain curve for resting muscle (schematic). 
Right: stress-strain curve for rubber (schematic). (Reproduced by permis- 
sion from original of figures 15a and 15b, E. Guth, Ann. NV. Y. Acad. Sci- 
ence, 47, 740). 


particle at r(z,y,z2) receives a displacement s(a,y,z) and the parti- 
cle at r+ dra displacement s+ ds, where ds = dr-Vs, then the 
relative position between the particles after the deformation is: 
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dr’ = dr+dr-Vs = dr-(I1 + D); 


I = idemfactor = S t,¢,» %, = unit vector along orthogonal 
= aXeS 2, Y, 2; (1) 


D=Vs; 


where Vs is the dyadic which gives the differential of strain at the 
point r(2,y, 2). 

Note that there is no restriction on the magnitude of s and hence 
none on Vs. The displacements and strain dyadics then may be of 
any finite magnitude. 


A 
at 
op 


wn 
S+dr-VS 
WN 


dr 


SS 


FicureE 5. Straining displacement of points in elastic medium during 
deformation. 


The dyadic D in equation (1) may be written in the usual way as 
the sum of a self conjugate and anti-self conjugate dyadic, D = 
S+Re= + (Vs + sV) + = (Vs—sV). Or in tensor notation: D,, = 
3(D;; + a) + = (D,,- D j;) =e,,+7?,;. Many writers on the sub- 
ject are not clear on the nature of the operation dr-D., In general, 
the argument is that in virtue of the decomposition D =§ + Ff, dr-S 
represents a ‘‘pure’’ strain and dr-F a rigid body rotation about an 
axis through the point P located at r(v,y,z2), referred to the unde- 
formed body. Thus, we find this argument for rigid body rotation by 
R in the usual works on elasticity, physics, and linear vector func- 
tions (Joos, 1932; Sommerfeld, 1950, p. 63; Wang, 1953; Weather- 
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burn, 1951, p. 138; etc.), and we even find mathematical reasons to 
support it (Gordon, 1950), albeit based on erroneous physical con- 
cepts. On the other hand, the nature of F as a straining rotation, 
as distinct from a rigid body rotation, is recognized by Swainger 
(1950) for axes in the deformed body and by Morse & Feshbach 
(1953, p. 67) for axes in the undeformed body. The physical bases 
given by the latter authors for the effect produced by F is readily 
recognized when it is explicitly stated that rigid body rotation can 
occur with reference to only one point of the undeformed body. 
Thus, any point P suffers a rigid rotation with reference to the 
fixed origin of axes 0. This is, however, not true of the point P’, 
separated by the vector distance dr from the point P when the latter 
is considered as origin. The meaning of dr-F is then clearly that 
of a straining rotation. When a medium is said to be deformed (the 
various points receiving relative displacements with respect to 
each other as origins), a rigid rotation can no longer obtain. Any 
rotation which occurs as part of the deformation must be in the 
nature of a straining rotation and not a rigid body rotation about the 
point P. The only alternative found in the strain produced by D is 
that of pure rotation unaccompanied by straining rotation, i.e., 
R=0, so that D=S. Since R may be written as -§ curl s x /, this 
necessitates the condition curl s =0 or the corollary that the dis- 
placement vector s may be written as a displacement potential 
s = Vd. The importance of recognizing FP as a definite contribution 
to the strain is obvious. 

We proceed now to give a formal proof to show that R cannot be 
a rigid body rotation and must of necessity be a straining rotation 
by deriving the necessary condition for rigid body rotation about 
the point P and showing immediately that the operation dr-R does 
not satisfy this requirement. Since in general FP ¢ 0, it follows that 
dr. R =-4 dr.(curl s x /) is a straining rotation. 

Proof. Let the elemental vector dr between two points fr and 
r+ dr, referred to orthogonal axes (2,y,2) in the undeformed body, 
be rotated in rigid fashion about some axis through P. Such an op- 
eration must leave the magnitude of dr unchanged. Let us repre- 
sent the operation on the vector dr by the specialized dyadic U 
with cartesian coordinates U; 3 %j = @,y,z. In order that the trans- 
formed vector dr. U should have the same magnitude as dr for any 
dr, we must have: 
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(dr. U)-(dr-U) = Je » Wann “ai 2h] 


m Fore! 


= ss pe u., “ni dr dr, = dr- dr = ve (dr)? , 


In m n 


1,m,n = 2,Y,2 


ljl=n 
2 ie Ee (2) 


m 


or 


In other words, the orthogonal relations for the components of the 
dyadic U must be satisfied for rigid body rotation. 
Consequently for R to represent rigid body rotation, it is neces- 


sary to satisfy = rit, = 5,,- Since R = 3(Vs - sV); Vs = 


| 
ds : oe Os | 
Da t,—— 3%, = %Y,2; 1, =4,),k; sV =)? eid 5 
an Ox 
n n # 


and, it follows at once that the orthogonal relations are not satis- 
fied. Since in general FR 40, the operation dr-R effects a change 
in both magnitude and direction of two deformable points relative 
to each other and is, therefore, a straining rotation of the point P’ 
relative to the point P. 

We know, moreover, from tensor analysis that the form of the 
dyadic D=8+R referred to the principal axes, when F& =0, is 


* ° . ° . ° 9 e 
,€,2, + tggt. + 76,7, where the e,’s are the eigenvalues. That 


is, the symmetric dyadic S referred to principal axis is ss te tne 
n 
In terms of the rectangular axes 2,y,z, when these are the prin- 
cipal axes, we have from equation (1): 


da =(1+e,)dz, 
dy’ =(1+e,)dy, (3) 
dz’ =(1+e,)dz, 


where dz is a line element parallel to the z axis before the deforma- 
tion, and dx’ the element after deformation, and e,, ey» &, the 
principal extensions. The quantities (1+e,), &=2,y,2 are the 
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magnifications of the line elements, and the values e, then meas- 
ure the increases in length per unit length in the directions of the 
principal axes. That is, they are the ‘‘strains’’ in these direc- 
tions. As already noted there is no restriction on their magnitudes. 
From equation (3), for a rectangular finite volume of sides a, B,+¢; 
the new dimensions of the rectangular volume after deformation, 
when referred to the principal axes 2,¥Y,2, are: 


a’ / (l+e,)dw, 
0 
b 


u (1 +e,)dy, (4) 


0 
¢ / (l+e,)dz. 
0 


In general, the quantities e, are functions of the dyadic compo- 
nents ¢,, and, hence, of the coordinates (#,y,z). When they are 
constant, equation (4) yields for the new dimensions: 


b’ 


oy 
I 


a=(1+e )a, 
b’ =(1+ e,)6 ; (5) 
c=(1+e,)c. 


If A is the cross-sectional area of the rectangle before the deforma- 
tion, the area after the deformation is: 


A’ = be’ =(1+ €,)(1 + e,) bo=(l+e,)(1+e,)A. (6) 


The theorem follows at once that A’ = A, if and only if the princi- 
pal strains e, are small compared to unity. 


IV. Stress-Strain Relations. 


In some texts on elasticity (Planck, loc. cit.; Sommerfeld, loc. 
cit.), the stress-strain relations are assumed to be linear since the 
Strains are taken to be small, and consequently, a linear form is 
obtained by expanding the stress in terms of the strains in a Tay- 
lor’s expansion and retaining only the first-order terms in the ex- 
pansion. Thus for example, in the most general form, referred to 
cartesian axes #,y,2, we have six relations of the type: 


he iF oeikt C12% yy + Cig@,, + C14 xy ve C15 eye ies eee 
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where Tr; is the stress acting in the direction of the j axis ona 
surface normal to the 7 axis, with 36 elastic constants which in 
virtue of showing that principal axes of strain are identical with 
principal axes of stress can be reduced to two independent elastic 
constants in an isotropic body. 

Indeed, the stress-strain relations will be linear if the strains 
are small. Such a restriction in the relation is, however, unneces- 
sary. Hooke’s law may be rendered freely as ‘‘Extension is pro 
portional to force’’ (Wang, loc. cit.), or ‘‘ut tensio, sic vis’’ (Joos, 
loc. cit.). In most cases it is a matter of common experience that 
Hooke’s law is no longer valid for large deformations, and so the 
strains are considered, perforce, small. 

In the following, we shall examine elastic behavior under the 
fundamental assumption that Hooke’s law applies to all bodies in 
the range of reversible elasticity, where the original force (stress) 
is replaced by ‘‘true’’ stress, to be defined below, and the original 
extension (strain) by strains of arbitrary magnitude. 

Under this assumption and the usual ones for an isotropic medium 
in equilibrium, the dyadic form of the stress-strain relation of true 
stress to strain of unrestricted magnitude is: 


Ted S11 + 28 (7) 


where T is the stress dyadic, S the a straining displacement 
component of the dyadic D, given by §,, 3(D;; + D,,) = ¢,;, and / 
the idemfactor. The factor | S| is the Lean: of the dyadic S, and 


is equal to 8 The symbols A,p denote, as usual, the Lamé 


constants of the medium. They are used for the relation between 
true stress and strain as well as when the deformation is of small 
order. It will be seen below, when stated in terms of the modulus 
of elasticity E and Poisson’s ratio v, that the meaning of these 
constants is now more general. The constants, however, acquire 
their usual more restricted meaning when the deformation is small. 
In tensor form, the equation is written simply as: 


T;; = A065, baal ot 
a-|sl, 7 ees 


O;74j 


(8) 
ey = Sy 0 = fet 
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Following Swainger (1947a,b) but using a reference system of 
orthogonal axes in the undeformed body rather than in the deformed 
body, the quantity T,, above is defined as the ‘‘true’’ stress, i.e., 
the load divided by the actual cross-section at the point r(z,y, 2). 
This-is a further generalization of Hooke’s law and is in fact im- 
plied by it since equation (8) in its most generalized form is a 
tensor relationship involving point functions. The alteration in 
area at the point is then very much a part of the law and can only 
be ignored in certain limiting cases such as when the deformation 
is small. 


V. Further Generalization of Hooke’s Law. 


We consider now, for convenience, a muscle fiber whose approxi- 
mate shape is that of a rectangle with sides parallel to the (a,y, 2) 
axes and with cross-section A (Figure 6). No distinction is made 


FIGURE 6. Schematic of muscle fiber cross-section A undergoing sim- 
ple elongation in the x direction and contraction in the y and 2 directions 
due to a normal force F acting in the positive 2x direction. 


between the nucleus and cytoplasm with respect to the elastic 
properties. The geometry of the muscle fiber will not invalidate 
the argument. Let the fiber be fixed at one cross-section, and ap- 
ply a normal force F’, = F, with F (constant) at the free end. The 
normal applied stress is then F_./A=F/A which is a simple 
tension. 

This normal stress is then a constant throughout each cross- 
section of the fiber. Since the ensuing deformation is now con- 
sidered to be relatively large, the area A undergoes a deformation 


as well. For a stress F/A in the 2 direction only, the equations 
(8) reduce to: 
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F 
quae eaters) + Oe pe], 


XX vy xx xx 
Ue (ena 6, te.) tren = Tl. , : (9) 
O=Al(e, + Pee é.,)t+2pe_. = Le, 


where A, are the elastic constants of the medium. The equations 
(9) show that the strain tensor for this case is 


Ort 0 0 

0 e 0 
vy 

0 0 e 


and is symmetric only. Thus, D = S + R is for this case D = S, i.e., 
R =0. Consequently, there is no rotational strain under these load 
conditions, as is to be expected, and moreover, the form of the 
stress and strain tensor show they are along the principal axes of 
the symmetric dyadic D=S. Thus e,, =e,, i al Pe Pes 
On physical grounds this is, of course, to be expected in the sim- 
ple tensile experiment since a normal force along one of the prin- 
cipal axes results only in pure strain introducing no shears or 
strain rotation. Consequently, from equation (6), the true stress is: 
F 

A(1+ e,) (1 + €,) 
to read: 


=€@ . 
z 


, and the first of equations (9) must be modified 


F 
A(1i + e,)(1 +e.) 
Equation (10) is the exact statement of Hooke’s law for elastic 
F 
Allie -e7 7) (1+ ea) ; : 
true stress at the point, equation (10) is a linear relation. The ap- 
proximate statement when the strains are small compared with unity 
is given in the first of equations (9) and is the limiting case usu- 
ally assumed in the analysis of engineering substances. 
From the second and third of equations (9), we find that: 
r 
=@€  =—-———eé_. 
YY ze 2(r re ut) X%% 


=A(e,+e,+¢,)+2pe,. (10) 


deformations. Writing = T*., where T’ , is the 


(11) 


A iret 
The ratio rAsepy of the transverse contraction to the longitudinal 
+ yl 
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elongation is Poisson’s ratio and is designated by v. The rela- 
tions between the Lamé constants A,p and the engineering con- 
stants E,v where E is Young’s modulus and v Poisson’s ratio are: 
(Sommerfeld) 


(3A + 2y) nN 
ee a oe YS eee 
Na? 2(A + p) 


= (12) 
Thus, using these definitions and equation (11), equation (10) may 
be written as: 
F 
aaa): ; (13) 


where we have dropped the subscript z, to indicate extension in the 
direction of the w axis. In terms of the true stress, equation (13) 
may be written as: 


TEE e; (14) 


Since the experiments in the literature, such as Ramsey and Street’s, 
represent the results of stress-strain relationships in terms of an 
initial cross-section area A (even when an average value is used) 
and not on a changing area depending on the actual deformation, 
the equivalent analytic expression for this representation is that 
given by expression (13). Hence, the result of experimental data 
on stress-strain diagrams when referred to the load per unit initial 
(static) area must then @ priori give the shape of a curve given by: 


F 
q Ee 2Ebve+ Eve’, (15) 


which is not a linear curve although a completely Hookian analysis 
has been used in its derivation. The only conditions used in ob- 
taining it are those of unrestricted strains and the consideration of 
the true stress when such is the case. Equation (14) defines 
Young’s modulus as T’/e = F/Ae(1- ve)’, i.e., the ‘‘true’’ stress 
divided by the strain. On the other hand, equation (15) allows the 
calculation of the constants E,v from the nonlinear curve. 

Equation (15) as it has been derived is a necessary condition for 
a generalized Hooke’s law involving true stress and strains of 
arbitrary magnitude. The proof that it is also a sufficient condition 
is easily given. In equation (15), the strains are of arbitrary mag- 
nitude and not ‘‘small,”’ since if the latter were true the equation 
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would immediately reduce to the classical equation F/A = Ee for 
small strains which is contrary to hypothesis. Writing equation (15) 
as F/A = Fe(1—ve)? and using the definition of Poisson’s ratio, 


— ve, =e, = e,,, equation (15) may be written as 


F/A(1 + e,,(l+e,,)=Ee or TT =Ee. 


Therefore, we have the result: 

The nonlinear stress-strain relation F/A = Ee(1-—ve)*, where F 
ts the applied normal force, A the initial area, e = normal strain of 
arbitrary magnitude, E = Young’s modulus (true stress divided by 
arbitrary strain), and v = Poisson’s ratio, is a necessary and suf- 
ficient condition for the generalized linear Hooke’s law. 
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FiIGuRE 7. Theoretical stress-strain curve of equation (15) for the 
range 0<e<1/v, v= 1/2, with stress referred to original area. Ordi- 
nate: force per unit initial area. Abscissa: increase in length per unit 
length. 


In the limiting case of ‘‘small’’ deformations, equation (15) re- 
duces to the well-known classical law of the simple tensile test, 
and Young’s modulus reduces to T/e = F/Ae as in the classical 
theory. The shape of the curve (15) is given in Figures 7, 8, and 


1 1 : 
9, for the ranges 0g e<—, @2-» and e>0 respectively. The 
v Vv 


164 WILLIAM H. GUTSTEIN 


curve of equation (15) has a maximum at (1/3 v, 4/27v) and a mini- 
mum (for which F/A = 0) at (1/v, 0). An inflection point is ob- 
served between the maximum and minimum at (2/3 v, 2/27v). The 
exact positions of these points depend on the magnitude of the 
modulus of linear elasticity and on Poisson’s ratio. Equation (15) 
is dimensionally correct because, since both e and v are pure num- 
bers and E is given in dynes per square centimeter in the c.g.s. 
system, both sides agree. From equation (12), since both A and p 
are always positive, v < + We shall assume the elastic properties 
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FIGURE 8. Theoretical stress-strain curve of equation (15) for the 
range ¢€ > 1/v, v= 1/2, with stress referred to original area. Ordinate: 
force per unit inztial area. Abscissa: increase in length per unit length. 


of the muscle fiber to be approximated by those of rubber, for which 
v is near the limiting value (Planck, Joc. cit.). Hence, in Figure 7 
the expression (15) is plotted with v = 5, for a strain range of 200% 
for purposes of comparison with the experimental curve of Ramsey 
and Street in the same range (Figure 1). In Figure 8 the same ex- 
pression is plotted for a strain range greater than 200% for pur- 
poses of comparison with Figure 2. In Figure 9 it is plotted for 


MUSCLE ELASTICITY 165 


the entire range of strain (juxtaposition of Figures 7 and 8) for com- 
parison with the total range curve of Ramsey and Street (Figure 3). 
With respect to the last curve (Figure 9), the range of e > 2 is ex- 
tended only for a relatively short distance to facilitate comparison 
with Figure 3, It may be noted that the experimental data of Ramsey 
and Street on active muscle (Figure 1) may be represented by a 
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FIGURE 9. Theoretical stress-strain curve of equation (15) for the 
range e > 0, v= 1/2, with stress referred to original area. Ordinate: force 
per unit inttial area. Abscissa: increase in length per unit length. 


curve having an inflection point near the upper end of the strain 
axis, similar to the theoretical curve of Figure 7. The solid line 
of Figure 1 is a curve given by an equation in Ramsey (Joc. cit.). 
Consequently, it may be seen that the experimental findings of 
Ramsey and Street on isolated muscle fibers may be accounted for 
in a qualitative fashion by a Hookian analysis, and even in a semi- 
quantitative fashion if the Poisson value v = 4 is assumed to hold 
for muscle. 

In the curves of Ramsey & Street, the resting length is defined 
as that length at which the fiber is found to develop maximal ten- 
sion upon stimulation. Reference to the theoretical argument shows 
that this corresponds to a strain value of e=1/3v. Thus a fiber 


1 cm in length develops a maximal tension at length cms. 
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However, according to the interpretation given by Ramsey & Street, 
the reason for this maximal developed tension is that the fiber at 
this length is under zero resting or initial tension. If the theoreti- 
cal curve, equation (15) is assumed to hold, then it becomes clear 
why the tension reaches a maximum at all. This is not due to zero 
resting tension at this length, but rather the natural result of relat- 
ing the true stress-strain curve to tension in terms of initial area. 
When the muscle fiber is passively stretched, it gives rise to an 
experimental curve of the type shown in Figure 2. Such a curve is 
again a natural result of equation (15) after the strain has reached 
a value of 1/v. An observation the meaning of which at first glance 
is not clear is that of the difference in the positions of the origins 
of this type of curve on the strain axes. This discrepancy can be 
explained in part as will be seen below. 

An interesting observation that can be made from Figure 3 of the 
total tension curve is that, after the maximal tension is reached, 
the ‘‘resting length’’ is exceeded and the effect of the rising pas- 
sive tension is brought to bear on the declining active (stimulated) 
tension of Figure 1. Despite this effect it does not impede the 
momentum of the latter chiefly because, as Figure 3 shows, it is of 
relatively small magnitude at first. However, after a strain value 
of 0.8 is reached the effects are noted in the upward sweep of the 
total tension curve. During the decline, total tension contains 
(algebraically) the two components of active and resting tensions. 
Yet from the standpoint of pure elastic phenomena, which concerns 
itself only with the relationship between tension and strain regard- 
less of individual contributions to either, the final result of the 
total tension curve is contained in the theoretical curve except for 
some modification in its shape and origin of axes. 

The conclusion must then obtain that equation (15) gives the 
proper form of the curve, which in its ascending portion corresponds 
to the experimental curve obtained by stimulation of the fiber alone, 
in its descending portion to the tension developed in the fiber by 
stimulation plus stretching, and finally in its upward sweep again 
to the natural development of a generalized stress-strain curve re- 
ferred to initial dimensions of the fiber. Equation (15) is simply 
the proper curve for the relationship between tension and strain. 
It does not discriminate between contributions to the tension. This 
information is provided by the separate experiments of Ramsey & 
Street. It does, however, tell why the curve has this peculiar 
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Shape. The latter derives from the reference system of stress. No 
attempt can be made to fit equation (15) to the experimental curves 
and to evaluate the constants E,v since, in principle, equation (15) 
demands that the original cross-sectional area A of the fiber be 
uniform. 

In experiment this is not the case since the area of the fiber is 
sometimes computed for a circular cylinder and at other times as 
an ellipse. At still other times average areas are used because the 
cross-sectional levels vary. In effect, this variation in the actual 
area may account for the failure of equation (15) to represent the 
experimental curves exactly. Of course, we have assumed that the 
muscle fiber behaves like an isotropic body throughout its range of 
elastic extension. Thus, it is not necessary to invoke the crystal- 
lization hypothesis to explain the peculiar shape of the stress- 
strain curve. This conclusion was also reached by Guth (loc. cit.). 


VI. Concluding Remarks. 


It is now clear that, although the evact xpression for an elastic 
substance of the relationship between stress and strain is the gen- 
eralized law of Hooke as long as it is within its elastic limit, this 
law will not have a linear appearance if the stress is expressed in 
terms of force per unit initial area in the case of highly deformable 
substances. This latter method of expression is only an approxi- 
mation of Hooke’s law and, as the analysis above shows, departs 
considerably from linearity in the case of large displacements. 
Considering the exact expression, however, we have either a linear 
or a nonlinear relation, depending on whether we use equation (14) 
or (15). In the exact expression, the true stress 7” is plotted 
against the strain (Figure 10). For this true stress, the slope of 
the curve is given by: dT’/de = E. 

The experimental curves found in the literature, referring as they 
do to highly deformable biological substances, offer poor approxi- 
mations to the ‘‘real’’ situations. Since this is the case, interpre- 
tations of these ‘‘complex’’ curves may offer problems which in 
reality do not exist. For example, the first component of the ap- 
proximate curve, Figure 1 or 7, shows a rise to a maximum in ten- 
sion with a subsequent decrease to zero tension, a type of behav- 
ior which is difficult to explain on physiologic grounds. However, 
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dT/de=E 


TRUE STRESS (T’) 


STRAIN (e) 


FIGURE 10. True stress-strain curve of equation (14). Ordinate: true 
stress. Abscissa: increase in length per unit length. 


this is an ‘‘apparent’’ situation.* The ‘‘real’’ situation is that the 
fiber does not conform to this picture at all since its stretch is not 
being referred to the true stress.’ The real situation is instead one 
of relative simplicity. The fiber is obeying a linear relation gen- 
eralized as in the discussion above. The same holds true for the 
second component of the curve, Figure 2 or 8, which again exhibits 
an apparently nonlinear ascent. 

In the preceding analysis we have attempted to show that if the 
primary assumption is made that the stress is everywhere propor- 
tional to the strain (Hooke’s law), where the magnitude of the lat- 
ter is of no consequence, as long as the substance is in its elastic 
limit, a nonlinear relationship must follow if stress-strain interpre- 
tation is made according to force per unit initial area. If the muscle 
fiber is assumed to behave elastically like rubber (v = 1/2), there is 
relatively good agreement with experimental results. 

If the generalization of Hooke’s law as has been done above im- 
plies nonlinearity (because of the failure to account for ‘‘true’’ 
stress), one cannot then use the argument of nonlinearity to negate 


*That is to say, viewed from theoretical grounds, the plot of strain 
against force per unit inztial area is a fiction or an artifice which does 
not reflect the real linear character of the relationship observed when the 
true stress is plotted against the strain. Since the theoretical analysis 
is an attempt to explain the experimental findings, the peculiar character 
of the experimental curve must be held to be nonreal in itself—the result 
of the special selection of coordinates chosen to express it. 
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Hooke’s law, unless ‘‘true’’ stress has been accounted for. This 
is, of course, a logical impossibility. That Hooke’s law may be 
assumed to hold for certain substances of large deformability has 
already been discussed above. It depends chiefly on the nature of 
the substance and confinement of the range of stretch to that of re- 
versible elasticity. 

Such a treatment and explanation of the phenomena involved does 
not attempt to contradict other treatments. What it does attemptisa 
necessity for inquiring into the nature of the ‘‘true’’ stresses under 
large deformations. ‘‘True’’ stresses cannot be the same as force 
per unit initial area, even if one does not assume the linear law to 
hold in the very beginning. Moreover, the treatment presented here 
is in the nature of a phenomenological analysis resting on a more 
proper mathematical development rather than a theory proceeding 
from a physical model to certain consequents. Such an attempt, 
for example, has been made by the use of elastomere theory. Thus 
the analysis given here does not and cannot come into conflict with 
any other theory which explains the facts, regardless of which ap- 
proximation is better. It can only be viewed as an alternative ex- 
position. It should, however, cause reflection about the concept of 
stress used in any theory. 

Viewed in terms of a generalized linear theory, a great simplifi- 
cation is offered to the cell biologist who must always be in query 
of the forces operating in various parts of the cell under normal and 
abnormal conditions. Since the elastic properties of protoplasm 
are among its most important ones (Seifriz, 1944), the applicability 
of a simplified theory to investigation of cell phenomena may prove 
fruitful. 


Summary. 

1. Experimental tension-length (stress-strain) curves on muscle 
fibers are reviewed with reference to their a posteriori implication 
concerning the failure of Hooke’s law to hold. 

2. Consideration of the true stress as a tensor point function in 
connection with the classical laws of elasticity allows for a further 
generalization of Hooke’s law. 

3. This generalization when restated in terms of the coordinates, 
which are equivalent to those used in the experimental data, re- 
sults in curves which predict the data in a qualitative and semi- 


quantitative manner. 
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4, The treatment is phenomenological with respect to elastic 
properties of muscle fibers and developed without reference to 
thermo-elastic or myophysiologic properties. 
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RECEIVED 12-1-55 


ANNOUNCEMENT OF NEW EXAMINATIONS 


THE U. S. CIVIL SERVICE COMMISSION, WASHINGTON, D. C. 


There is an urgent need for Chemists, Mathematicians, Metallur- 
gists, Physicists, and Electronic Scientists in the Washington, 
D. C. area, the United States Civil Service Commission has an- 
nounced. Vacancies are in various Federal agencies and pay sal- 
aries ranging from $4,345 to $11,610 a year. 

To qualify for positions paying $4,345 a year, applicants must 
have had appropriate education or a combination of education and 
experience. For the position of Electronic Scientist, appropriate 
technical or scientific experience alone may be qualifying. For 
higher-grade positions, professional experience is also required. 
Graduate study may be substituted for all or part of this experience, 
depending on the grade of position. No written test is required. 

Further information and application forms may be obtained at 
many post offices throughout the country, or by writing to the U. S. 
Civil Service Commission, Washington 25, D. C. Applicants should 
ask for Announcement No. 46(B). Applications will be accepted by 
the Board of U. S. Civil Service Examiners, National Bureau of 
Standards, Washington 25, D. C., until further notice. 
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